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Chief,  Technical  Inf oraation  Division 

This  Final  Report  covers  all  work  performed  under  grant  number  AFOSR* 
82-0079,  covering  the  period  15  January  1982-14  January  1985.  The  overall 
goal  of  this  study  was  to  identify  and  evaluate  the  environmental  or 
microphysical  parameters  that  control  the  efficiency  of  the  various  mechanisms 
of  precipitation  development.  Such  evaluations  can  then  be  used  as  the  basis 
for  studying,  or  even  predicting,  the  effect  of  geographical  or  climatological 
differences  between  regions  on  the  microphysical  mechanisms  of  precipitation 
development.  f 

This  study  includes  work  on  warm  rain  initiation  and  development, 
ice  multiplication,  snowflake  aggregation,  and  the  growth  of  graupel  by 
collection  of  supersooled  water  droplets.  During  the  first  year  of  the 
study,  work  concentrated  on  the  warm-cloud  studies  and  on  ice  multiplication. 
During  the  second  year,  work  continued  on  the  warm-cloud  studies  and  began 
on  the  investigation  of  snowflake  aggregation.  With  the  extension  of  the 
study  into  a  third  year,  additional  investigations  were  made  on  showflake 
aggregation  and  preliminary  investigations  were  begun  to  extend  these  studies 
into  mixed-phase  (riming)  precipitation.  The  studies  of  riming  and  graupel 
development  were  quite  successful  and  indicated  that  the  studies  of  warm 


rain  initiation  and  development  could  be  extended  to  include  mixed-phase 
precipitation  development.  These  background  studies  have  been  summarized 
in  a  article  being  submitted  for  publication,  which  is  presented  in  the 
appendix.  Other  articles  resulting  from  this  study  are  under  preparation 


and  should  be  completed  in  the  near  future. 

In  general,  the  study  identifies  a  number  of  key  parameters  that 
control  the  microphysical  development  of  precipitation.  The  two  primary 
parameters  that  need  to  be  measured  on  a  climatological  basis  are  cloud  base 
temperature  and  cloud  droplet  concentration.  With  a  knowledge  of  these 
parameters  it  is  possible  to  estimate  the  liquid  water  content  and  mean 
droplet  diameter  as  functions  of  the  depth  of  a  cloud.  This  in  turn  permits 
an  estimate  of  the  likelyhood  of  warm  rain  development,  ice  multiplication, 
and  ice  crystal  concentration  in  the  cloud  (see  Hobbs  and  Rangno,  1985: 
Journal  of  the  Atmospheric  Sciences ,  42 ,  2523-2549. ).  Combined  with 
satellite  observed  cloud  tops,  a  climatological  knowledge  of  cloud  base 
temperature  and  droplet  concentration  would  go  a  long  way  toward  predicting 
the  type  of  particle  development  likely  in  widely  different  geographical 
areas.  Other  key  parameters  include  the  degree  of  entrainment  and  stability 
of  the  environment.  Entrainment  is  still  difficult  to  predict  in  any  quant¬ 
itative  way,  but  may  not  be  as  variable  geographically  as  some  of  the  other 
parameters.  Studies  of  stability  are  currently  possible  using  the  large 
world-wide  data  base  of  sounding  data  currently  available.  Cloud  base 
temperature  and  droplet  concentration  data,  on  the  other  hand,  are  still 
rare  for  most  parts  of  the  world. 
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BACKGROUND  STUDIES 


From  the  very  onset  of  these  studies  it  has  been  clear  that  cloud 
base  temperature  is  a  uniquely  important  parameter  for  assessing  geographical 
variability  in  the  microphysical  mechanisms  of  precipitation  development.  .It 
♦  has  significance  in  areas  ranging  from  the  overall  water  budget  and 
precipitation  efficiency  to  the  actual  details  of  the  cloud  microstructure 
and  its  evolution.  If  entrainment  is  neglected,  for  example,  then  the 
consensationally  produced  liquid  water  contents  at  any  height  above  cloud 
base  are  uniquely  defined  by  the  temperature  and  pressure  at  cloud  base. 

The  build-up  of  condensed  water  with  height,  of  course,  is  the  basic  driving 
force  behind  precipitation  development.  Estimates  of  these  adiabatic  water 
contents  are  shown  in  Figure  1  (after  Ludlam,  1950).  In  this  presentation, 
the  cloud  base  pressures  have  been  varied  in  parallel  with  the  cloud  base 
temperatures  in  order  to  produce  a  single  graph  having  applicability  to  a 
wide  variety  of  summertime  convective  clouds.  The  water  content  profiles, 
however,  are  not  very  sensitive  to  variations  in  cloud  base  pressure.  This 
means  that  a  single  figure  of  this  type  can  be  used  to  estimate  water  contents 
for  many  different  sets  of  meteorological  conditions. 

In  a  sense.  Figure  1  represents  the  maximum  liquid  water  content  that 
would  be  expected  in  the  absence  of  a  precipitation-related  accumulation 
zone.  The  average  water  contents  that  are  observed  in  convective  clouds, 
however,  are  typically  much  lower  than  the  adiabatic  values  illustrated  in 
Figure  1.  Warner~(.1970},”  for  example,  summarized  a~nuraber  of  observations 

of  the  ratio  of  the  observed  mean  liquid  water  content  at  a  given  height 
to  the  corresponding  adiabatic  value.  Although  there  is  considerable  variation 
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Figure  1.  Adiabatic  liquid  water  contents  for  cumuliform  clouds. 
Saturation  with  respect  to  liquid  water  is  assumed 
throughout  with  the  isopleths  iabeled  in  g  m"3.  The 
dashed  lines  indicate  the  corresponding  temperatures 
within  the  clouds. 


in  the  observations  there  is  a  consistent  trend  with  the  ratio  decreasing 
with  increasing  height  above  cloud  base.  This  decrease  can  be  approximated 
by  an  empirical  relation 


F  =  Q.2  + 


0.16 
H  +  .20 


(1) 


where  H  is  the  height  above  cloud  base  in  km  and 


L  =  F  La  ,  (2) 

where  L  and  L&  are  the  mean  liquid  water  content  and  adiabatic  liquid  water 
contents  respectively.  Figure  2  illustrates  the  effect  of  applying  (1)  to 
the  water  contents  shown  in  Figure  1 . to  produce  an  estimate  of  the  expected 
mean  water  contents  as  a  function  of  height  above  cloud  base. 

Although  cloud  base  temperature  has  been  frequently  ignored  in 
summarizing  project  results,  I  was  able  to  compile  some  data  on  base 
temperatures  of  summertime  convective  clouds  in  the  United  States  and  Canada. 
Figure  3a  shows  the  observed  distributions  of  cloud  base  temperatures  for 
Florida  and  Hawaii,  while  3b  shows  the  distributions  for  the  Midwest  and  the 
Great  Plains.  The  Florida  data  was  obtained  during  the  FACE-2  experiment  in 
1978,  1979,  and  1980  and  includes,  data  from  48  different  operational  days. 
Hawaiian  data  were  obtained  from  Professor  Takahashi  at  the  Cloud  Physics 
Observatory  in  Hilo  and  represent  observations  taken  over  32  different  days. 
The  Midwest  data  was  combined  from  three  different  projects  to  give  data  for 
72  days.  The  Great  Plains  data  was  similarly  combined  from  the  NHRE,  HIPLEX, 
and  CCOPE  projects  to  give  data  for  123  different  days.  Recently,  similar 
data  was  obtained  from  Dr.  Schemenauer  at  the  Atmospheric  Environment  Service 
(Canada)  for  two  additional  Canadian  sites  (Fig.  3c) . 
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Figure  2.  Average  liquid  water  contents  for  cuniuliform  clouds. 

Saturation  with  respect  to  liquid  water  is  assumed 
throughout  with  the  isopleths  labeled  in  g  The 

dashed  lines  indicate  the  corresponding  temperatures 
within  the  clouds. 
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Figure  3a.  Observed  distributions  of  cloud  base  temperature 
for  summertime  convective  clouds  in  Hawaii  and 
Florida. 
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Figure  3b.  Observed  distributions  of  cloud  base  temperature  for 
summertime  convective  clouds  in  the  Great  Plains  and 
the  Midwest. 
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Figure  3c.  Observed  distributions  of  cloud  base  temperature  for 
summertime  convective  clouds  in  Montana  (Miles  City) 
and  two  Canadian  sites  (Thunder  Bay  and  Yellowknife) 


These  distributions  are  rather  interesting  in  that  they  highlight  the 
wide  variety  of  cloud  base  temperatures  expected  in  summertime  convective 
clouds.  Perhaps  the  most  interesting  feature  of  this  data  is  how  moist 
Midwest  clouds  can  be,  often. having  temperatures  in  the  same  range  as  found 
in  Florida  or  Hawaii.  Cloud  base  temperatures  in  the  Great  Plains,  on  tha 
other  hand,  are  significantly  colder  than  those  in  the  other  two  regions 
with  little  overlap  between  the  observed  distributions.  Studies  of  clouds 
in  the  Great  Plains  are  thus  not  likely  to  be  relevant  to  clouds  over  much 

« 

of  Central  and  Eastern  areas  of  the  United  States. 

WARM  RAIN  INITIATION 

In  earlier  studies  (e-g-»  Johnson,  1978;  1980)  it  was  suggested  that 
a  quantitative  threshold  could  be  defined  that  would  signal  the  onset  of 
effective  coalescence  growth.  This  concept  was  the  direct  result  of  a  - 
limited  number  of  tests,  with  the  threshold  that  was  ultimately  suggested 
^  >  4  x  10-3  sec-l)  growing  directly  out  of  a  graphical  examination 
of  drop  growth  as  a  function  of  time.  In  the  current  studies,  this  concept 
was  re-examined  and  tested  over  a  wide  range  of  conditions.  In  these 
studies,  possible  thresholds  were  redefined  in  terms  of  the  mass  doubling 
time  of  the  growing  drop.  Two  different  thresholds  were  examined  -  mass 
doubling  in  3  minutes  and  mass  doubling  in  5  minutes.  The  mass  doubling 
time  (MOT)  is  a  more  appealing  way  of  expressing  the  coalescence  threshold 
since  it  can  be  relatively  easily  interpreted  in  terms  of  significant  drop 
growth  over  cloud  lifetimes.  Since  large  cloud  droplets  beginning  their 
growth  to  raindrops  will  have  to  double  in  mass  approximately  ten  times 
before  reaching  raindrop  proportions,  these  possible  thresholds  effectively 


restrict  the  description  of  growth,  as  "effective"  to  those  situations  in 
which  coalescence  rain  could  be  produced  in  the  30-50  min  lifetime  of  a 
convective  cloud. 

Figures  4-11  illustrate  the  growth  of  large  drops  in  a  variety  of 
situations.  In  each  case,  drop  growth  is  modeled  in  a  constant  updraft 
framework  in  which  continued  condensation  results  in  a  steady  build  up 
in  water  content  as  the  parcel  is  lifted.  In  a  separate  series  of  runs 
drop  growth  was  computed  in  clouds  having  a  constant  liquid  water  content 
and  negligible  updraft  (1  cm  s”*)  .  In  effect,  the  former  runs  correspond 
to  growth  in  a  convective  cloud,  while  the  latter  runs  simulate  stratiform 
clouds. 

The  model  used  for  these  calculations  is  essentially  the  continuous 
collection  parcel  model  discussed  in  Johnson  (1982).  For  use  in  this 
study,  however,  the  model,  was  modified  ot  include  an  adjustment  in  the 
calculated  terminal  velocities  of  cloud  droplets  and  raindrops  to  account 
for  changes  in  air  density  as  the  parcel  rises.  In  addition,  the  model 
was  modified  to  permit  drop  growth  by  condensation  to  be  arbitrarily 
restricted  to  any  fraction  of  the  expected  growth  in  an  adiabatic  parcel. 

If  desired,  the  condensational  growth  can  also  be  controlled  to  yield  water 
contents  that  follow  the  empirically  derived  mean  water  contents  specified 
by  equation  (1) . 

In  the  convective  cases  (Figs.  4-11)  the  basic  pattern  is  one  of 

slow  initial  growth  followed  at  some  point  by  a  sudden  spurt  of  growth. 

This  overall  pattern  and  the  seeming  suddenness  of  the  transition  from 

«* 

slow  growth  to  fast  growth  were  the  original  features  which  suggested  that 
it  might  be  possible  to  define  the  transition  from  one  growth  regime  to 
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Figure  6.  Growth  of  large  drops  (30-200  ym  initial  diameter)  in  a  simple 
cloud  model.  The  upper  plot  assumed  a  cloud  base  temperature 
of  0°C.  The  middle  plot  assumed  a  10°C  base  temperature  and 
the  bottom  plot  assumed  a  base  temperature  of  20°C.  Each  case 
involved  adiabatic  water  contents. 
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Figure  7.  Growth  of  large  drops  (30-200  yra  initial  diameter)  in  a  simple 
cloud  model.  The  upper  plot  assumed  adiabatic  water  contents, 
the  middle  plot  assumed  water  contents  that  were  half  of  the 
adiabatic  values,  and  the  bottom  plot  assumed  a  height  dependent 
fraction  of  the  adiabatic  water  content  that  corresponds  to 
Figure  2.  In  each  case,  the  cloud  base  temperature  was 
assumed  to  be  10°C. 
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Figure  8.  Growth  of  large  drops  (30-200  Jim  initial  diameter)  in  a  simple 
cloud  model.  The  upper  plot  assumed  a  cloud  base  temperature 
of  0°C.  The  middle  plot  assumed  a  10°C  base  temperature  and 
the  bottom  plot  assumed  a  base  temperature  of  20°C.  Each  case 
involved  adiabatic  water  contents. 
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Figure  9.  Growth  of  large  drops  (30-200  yra  initial  diameter)  in  a  simple 
cloud  model.  The  upper  plot  assumed  adiabatic  water  contents, 
the  middle  plot  assumed  water  contents  that  were  half  of  the 
adiabatic  values,  and  the  bottom  plot  assumed  a  height  dependent 
fraction  of  the  adiabatic  water  content  that  corresponds  to 
Figure  2.  In  each  case,  the  cloud  base  temperature  was 
assumed  to  be  10°C. 
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Figure  11.  Growth  of  large  drops  (30-200  yra  initial  diameter)  in  a  simple 
cloud  model.  The  tipper  plot  assumed  adiabatic  water  contents, 
the  middle  plot  assumed  water  contents  that  were  half  of  the 
adiabatic  values,  and  the  bottom  plot  assumed  a  height  dependent 
fraction  of  the  adiabatic  water  content  that  corresponds  to 
Figure  2.  In  each  case,  the  cloud  base  temperature  was 
assumed  to  be  10°C. 
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the  other  in  terms  of  a  "coalescence  threshold."  The  symbols  superimposed 
on  each  growth  curve  indicate  the  points  at  which  the  coalescence  growth 
alone  is  adequate  to  make  a  MDT  of  5  minutes  (left-most  symbol)  or  3 
minutes  (right-most  symbol) .  Either  of  these  possible  thresholds  is 
capable  of  partitioning  the  growth  curves  into  a  slow  initial  growth 
regime  and  a  rapid  growth,  regime  that  quickly  produces  drops  of  raindrop 
dimensions.  Of  these  two  possible  thresholds,  the  3-min  doubling  is  a 
rather  conservative  measure  of  the  onset  of  rapid  growth  while  the  5-min 
doubling  is  a  rather  optimistic  estimate  of  the  transition  point.  In 
most  cases,  however,  the  transition  from  slow  to  rapid  growth  is  so  fast 
that  these  two  alternatives. for  defining  a  threshold  value  are  not  very- 
far  apart.  This  is  particularly  true  in  the  cases  for  which  drop  growth 
is  relatively  rapid.  For  this  reason,  the  most  sensitive  comparisons  of 
the  two  possible  thresholds  are  those  runs  in  which  the  evolution  of 
raindrops  is  relatively  slow  -  i.e.,  those  runs  which  have  the  coolest 
cloud  base  temperatures,  highest  total  drop  concentrations,  weakest 
updrafts,  or  most  entrainment. 

The  stratiform  runs  examined  drop  growth  in  situations  in  which  the 
weak  updrafts  present  do  not  force  the  cloud  development  by  the  continual 
buildup  in  condensed  liquid  water,  but  rather  have  liquid  water  contents 
that  are  effectively  constant  with.  time.  In  these  cases  the  onset  of 
effective  coalescence  growth  was  less  dramatic  and  the  critical  parameter 
for  raindrop  or  drizzle  production  seemed  to  be  the  overall  time  available 
for  growth  to  continue.  Although  the  coalescence  threshold  (MDT  of  3  or  5 
minutes)  didn’t  signify  any  marked  change  in  the  stratiform  growth  curves, 
it  still  serves  as  a  rough  indicator  of  the  relative  speed  of  drop  growth. 


To  be  useful  in  the  context  of  interpreting  cloud  structure  and 
evolution,  the  coalescence  threshold  (MDT  of  3-5  minutes)  must  be 
re-expressed  in  terms  that  directly  relate  to  cloud  thickness  or  degree 
of  development.  Figures  12  and  13  are  a  start  in  this  direction, 
identifying  those  drops  that  will  grow  at  or  above  the  3  (Fig.  12)  or 
5  (Fig.  13)  minute  MDT  threshold  for  liquid  water  contents  between  0  and 
3  grams  per  cubic  meter.  In  each  Figure,  plots  have  been  generated  for 
each  of  six  different  total  droplet  concentrations  (N  =  50,  100,  200, 

400,  800,  and  1200  cnf  ^) .  It  is  clear  that  there  are  minimum  liquid 
water  contents  (and  minimum  drop  diameters)  for  effective  coalescence 
growth.  If  the  shaded  region  of  each  plot  were  strictly  rectangular,  it 
would  be  possible  to  identify  specific  minimum  water  contents  and  drop 
diameters  for  effective .coalescence  growth.  Since  the  lower- left -hand 
comer  of  the  shaded  region  is  curved,  however,  there  is  a  degree  of 
uncertainty  in  assigning  any  single  minimum  value  of  liquid  water  content 
as  marking  the  transition  to  effective  coalescence  growth.  One 
identifiable  critical . water  content  is  the  lowest  value  for  which  any 
size  drop  can  meet  the  threshold  criteria  of  a  3  or  5  minute  MDT.  In 
general,  however,  drops  meeting  this  criteria  will  have  to  already  be 
approaching  millimeter  dimensions .  The  more  crucial  question  is  when 
drops  in  the  vicinity  of  50  pm  diameter  are  capable  of  starting  effective 
coalescence  growth  since  these  drops  may  be  numerous  enough  to  have 
significant  impact  on  subsequent  cloud  development  (Johnson,  1982).  This 
means  that  some  compromise  must  be  invoked  between  the  minimum  water  content 
that  can  support  effective  coalescence  for  any  sized  drop  and  selection  of 
a  water  content  that  permits  the  smallest  diameter  (and  hence  the  most 
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Plots  identifying  the  liquid  water  contents  (shaded  region) 
which  are  capable  of  supporting  effective  coalescence  growth 
(MDT  =  3  minutes)  as  a  function  of  the  diameter  of  the 
growing  drop. 


DROP  DIAMETER  (fim) 


Plots  identifying  the  liquid  water  contents  (shaded  region) 
which  are  capable  of  supporting  effective  coalescence  growth 
(MDT  =  5  minutes)  as  a  function  of  the  diameter  of  the 
growing  drop. 


numerous)  drops  possible  to  grow  effectively.  Close  examination  of  the 
plots  included  in  Figures  12  and  13  (as  well  as  other  plots  not  shown)  led 
to  a  simple  relation  that  may  be  adequate  to  specify  a  rough  estimate  of 
the  minimum  water  content  that  could  support  effective  coalescence  growth 
as  a  function  of  the  total  droplet  concentration: 

Lmin  *  °'6  *  °-01  N°'S  '  (3) 

where  Lm^n  is  the  minimum  liquid  water  content  in  g  m-3  and  N  is  the  total 
droplet  concentration  in  cm*3.-  This  particular  relation  corresponds  rather 
well  to  both,  the  minimum  value  of  water  content  that  can  ever  meet  a  MDT 
threshold  of  5  minutes  and  (simultaneously)  allow. relatively  small  large 
drops  to  exceed  a  MDT  threshold  of  3  minutes. 

A  similar,  but  independent,  estimate  of  the  minimum  radius  of  large 
drop  capable  of  effective  coalescence  growth  led  to  the  relation 

logia  Rmin  *  15  *  °-006  N°’S  C«) 

where,  as  before,  N  is  the  total  droplet  concentration  in  cm*3  and  R  .  is 

mxn 

the  minimum  large  drop  radius  in  microns. 

Figure  14  repeats  several  of  the  growth  curves  for  individual  large 
drops  that  were  previously  discussed,  hut  with  additional  shading  to 
identify  the  region  of  the  plot  which  (3)  and  (4)  identify  as  being  capable 
of  supporting  effective  coalescence  growth.  While  the  correspondence  with 
the  individual  growth  curves  is  not  perfect,  it  is  apparent  that  these 
suggested  relations  can  be  useful  tools  in  quantifying  the  start  of  effective 
coalescence  growth. 

Another  way  of  looking  at  this  same  threshold  is  illustrated  in 
Figure  IS.  This  figure  is  a  copy  of  Figure  2,  but  with  an  overlay  identifying 
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15.  Mean  liquid  water  contents  for  cumuliform  clouds 
superimposed  with  shading  to  indicate  the  regions 
which  might  be  expected  to  support  effective 
coalescence  growth. 
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the  region  of  effective  coalescence  growth  as  defined  by  (3)  for  a  cloud 
base  droplet  concentration  of  600  cm-3.  This  figure  suggest  a  major  difference 
between  clouds  with  base  temperatures  substantially  below  10°C  and  those  with 
base  temperatures  substantially  above  10°C.  In  either  case,  effective 
coalescence  can  begin  within  an  adiabatic  core  within  a  kilometer  or  so 
of  cloud  base.  In  the  warm  based  clouds,  however,  effective  coalescence 
growth  can  still  occur  even  in  regions  that  have  average  water  contents. 

On  the  other  hand,  in  cold-based  clouds  average  water  contents  would  not 
be  adequate  to  support  effective  coalescence  growth. 

TIME  AND  LOCATION  OF  WARM  RAIN  DEVELOPMENT 

Following  the  study  of  coalescence  onset,  additional  studies  were 
performed  to  identify  the  time  and  location  of  precipitation  formation  as 
a  function  of  the  cloud  base  temperature,  updraft  velocity,  total  droplet 
concentration,  and  degree  of  entrainment.  These  studies  made  up  the  bulk 
of  the  first  year  effort  of  this  study.  Ihe  initial  studies  of  this  type 
all  used  the  same  continuous  collection  parcel  model  employed  in  the  studies 
of  the  growth  of  individual  particles  (Figs.  4-11). 

Table  1  shows  some  of  the  results  of  these  studies  -  presenting  the 
height  (m)  at  which  the  radars  reflectivity  factor  first  exceeds  10  dBZ. 

This  is  a  commonly  used  measure  of  precipitation  onset  and  corresponds  to 
the  first  development  of  significant  concentrations  of  raindrops.  In  all 
cases  in  which  the  parcel  exceeded  8  km  altitude  or  fell  below  -40°C  before 
reaching  10  dBZ,  the  height  entry  is  replaced  with  a  full  field  of  asterisks. 
Tables  2  and  3  present  similar  data,  but  in  these  tables  sub-adiabatic  water 
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contents  were  employed  (half  of  adiabatic  in  Table  2  and  a  variable 
fraction  of  adiabatic  corresponding  to  equation  (1)  in  Table  3).  One 
of  the  most  interesting  results  of  this  study  is  the  relative  insensitivity 
of  the  time  and  location  of  precipitation  onset  to  total  droplet  concentration. 

The  first  echo  heights  can  be  redefined  in  terms  of  the  additional  time 
Cor  height)  required  to  reach  the  10  dBZ  echo  intensity  after  passing  the 
threshold  water  content  defined  by  (3).  These  "adjusted"  heights  are 
presented  in  Tables  4,  5,  and  6  for  the  same  three  bulk  entrainment  rates 
invoked  in  generating  Tables  1-3.  Adjustments  of  this  sort  to  allow  for 
the  initial  "incubation"  period  of  a  cloud  in  which  coalescence  growth  is 
either  absent  or  ineffective  may  ultimately  simplify  predictions  of  the  time 
and  location  of  precipitation  onset  and  reduce  some  of  the  variability 
introduced  by  differing  updraft  velocities  or  entrainment  rates. 

The  parcel  model  used  for  these  studies  is  an  attractive  tool  for 
investigations  of  this  sort.  It  is  relatively  inexpensive  to  run  and 
produces  a  rather  precise  estimate  of  the  height  of  precipitation  development. 
These  are  both  valuable  attributes  when  performing  a  large  set  of  runs  for 
intercomparison.  While  intercomparisons  between  runs  are  relatively  reliable 
with  this  model,  the  absolute  accuracy  of  the  locations  of  precipitation  onset 
is  suspect  since  the  model  neglects  differential  sedimentation  of  the  falling 
drops.  This  problem  can  be  addressed  by  using  a  "trajectory"  model  that 
allows  explicit  consideration  of  the  relative  motion  of  the  growing  raindrops. 
When  incorporated  into  an  appropriate  framework  (see  Johnson,  1982) ,  this 
sort  of  model  can  produce  time-height  cross-sections  of  radar  reflectivity 
that  behave  very  much  like  observed  radar  "first  echoes"  (e.g.,  Battan,  1953). 
Figure  16  shows  one  such  cross-section  calculated  with  an  improved  version 
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Figure  16.  Tine-height  reflectivity  cross-section  from 

"trajectory4'  model  of  precipitation  initiation. 
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of  the  "trajectory"  model.  Rather  than  producing  a  single  estimate  of 
the  height  at  which  the  radar  reflectivity  reaches  detectable  limits,  this 
sort  of  model  actually  predicts  the  base  and  top  heights  of  the  echb  as  a 
function  of  time.  By  scanning  over  the  calculated  cross-sections,  it  is 
then  possible  to  estimate  the  average  heights  of  the  echo  top  and  base 
expected  in  an  observational  study. 

The  current  version  of  this  model  has  been  improved  to  include  the 
velocity  adjustments,  for  changing  air  densities  and  variable  entrainment 
algorithms  already  introduced  into  the  parcel  model.  In  addition,  in 
order  to  improve  the  resolution  of  the  calculated  cross-sections,  it  has 
been  necessary  to  increase  the  number  of  drop  categories  used  in  the 
computations  by  a  factor  of  four  over  that  used  in  the  parcel  calculations 

(to  101  categories).  The  extension  of  the  warm  rain  studies  to  include  results 
from  the  trajectory  model  comprised  the  major  portion  of  the  year  two  efforts 
of  this  study. 

A  full  set  of  runs  of  this  model  have  been  completed,  almost  a  thousand 
in  all,  covering  a  wide  range  of  updraft  velocities,  cloud  base  temperatures  (and 
pressures),  droplet  concentrations,  and  degrees  of  entrainment.  Table  7  shows 
the  temperature  of  the  predicted  first  echo  for  an  adiabatic  updraft  (°C). 
Temperatures  for  both  the  top  and  the  base  of  the  10  dBZ  3-minute  first  echo 
are  shown.  [The  "three-minute  first  echo"  is  just  the  average  top  and  base  of 
the  echo  which  has  reached  a  10  dBZ  threshold,  followed  ever  three  minutes  of 
development.]  At  the  upper  left  hand  side  of  the  table,  no, data  are  given  since 
the  model  never  developed  a  10  dBZ  echo.  At  the  lower  left,  one  run  is  also 
indeteminant  since  no  echo  was  formed  within  one  hour ,  and  calculations  were 


TABLE  7 


10  DBZ  FIRST  ECHO  TEMPERATURE  (TOP/BASE) 


ADIBAT  -  1.0 


CLOUD  BASE  TEMPERATURE 


UPDRAFT  VELOCITY,  M/S 


(8  CONCENTRATION,  CM-3) 

-5 

0 

5 

10 

15 

20 

25 

30 

16. C 

(  50.) 

««*/*** 

»»*/»** 

-44/-30 

-36/-17 

-23/  -8 

-11/ 

1 

0/  8 

7/ 

16 

16.0 

(  200.) 

***/*** 

••*/*** 

-43/-30 

-33/-17 

-20/  -8 

-8/ 

1 

0/  9 

9/ 

16 

16.0 

(  600.) 

ft** /*•* 

***/*** 

-44/-32 

-36/-18 

-20/  -8 

-8/ 

1 

0/  8 

9/ 

15 

16.0 

(1200.) 

ft**/*** 

/*•* 

-42/-34 

-35/-20 

-23/  -9 

-11/ 

0 

0/  8 

7/ 

IS 

14.0 

(  50.) 

*** /*** 

-40/-38 

-40/-25 

-2B/-14 

-17/  -5 

-6/ 

3 

2/  11 

11/ 

17 

14.0 

(  200.) 

•**/••* 

-40/- 39 

-41/-24 

-29/-13 

-17/  -4 

-6/ 

4 

3/  11 

10/ 

17 

14.0 

(  600.) 

*** /**• 

***/*** 

-40/-27 

-28/-15 

-17/  -6 

-6/ 

3 

1/  10 

10/ 

17 

14.0 

(1200.) 

*••/*** 

/*** 

-40/-29 

—31/— 16 

-16/  -7 

-8/ 

3 

2/  10 

11/ 

16 

12.0 

(  50.) 

***/*** 

-39/-30 

-36/-19 

-24/-10 

-13/  -2 

-3/ 

6 

4/  12 

12/ 

18 

12.0 

(  200.) 

***/*•« 

-40/-32 

-34/-20 

-24/-10 

-14/  -1 

-3/ 

6 

5/  12 

12/ 

19 

12.0 

(  600.) 

•ft*/*** 

—39/— 33 

-36/-21 

-24/-U 

-13/  -2 

-3/ 

5 

6/  12 

12/ 

19 

12.0 

(1200.) 

•** /*** 

-39/-3S 

—37/— 2  3 

-26/-12 

-12/  -3 

-3/ 

4 

4/  12 

12/ 

18 

10.0 

(  50.) 

-38/-31 

-40/-23 

-31/-14 

-18/  -6 

-9/  1 

0/ 

8 

8/  14 

14/ 

21 

10.0 

(  200.) 

-3S/-33 

-40/-23 

-31/-14 

-19/  -6 

-9/  2 

0/ 

8 

8/  IS 

14/ 

21 

10.0 

(  600.) 

-36/-3S 

-39Z-26 

-30/-16 

-18/  -7 

-9/  1 

0/ 

8 

8/  14 

14/  21 

10.0 

(1200.) 

«**/**• 

-39/-2B 

-31/-17 

-20/  -8 

-11/  0 

-2/ 

7 

7/  14 

15/  20 

8.0 

(  50.) 

—40/— 2  7 

-35/-20 

-25/-10 

-14/  -2 

-5/  4 

3/ 

11 

10/  17 

16/  23 

8.0 

(  200.) 

-38/-2S 

-36/-17 

-25/  -9 

-14/  -2 

-4/  4 

4/ 

10 

10/  17 

17/ 

22 

8.0 

(  600.) 

-38/-29 

-37/-19 

-24/-12 

-13/  -4 

-5/  4 

3/ 

10 

10/  16 

16/ 

22 

8.0 

(1200.) 

-38/-31 

—36/— 2 1 

-27/-12 

-14/  -4 

-7/  3 

1/ 

10 

9/  16 

16/ 

22 

6.0 

(  50.) 

—  37/— 2  2 

-28/-14 

-17/  -7 

-8/  0 

-1/  7 

6/ 

13 

13/  19 

19/ 

25 

6.0 

(  200.) 

-38/-21 

-27/-14 

-17/  -7 

-8/  0 

-1/  7 

6/ 

13 

13/  19 

20/ 

24 

6.0 

(  600.) 

-37/-20 

-29/-14 

-20/  -6 

-9/  0 

-2/  7 

6/ 

12 

12/  19 

19/ 

25 

6.0 

(1200.) 

-37/-20 

-30/-14 

-19/  -7 

-11/  1 

-3/  7 

6/ 

12 

13/  18 

18/ 

24 

4.0 

(  50.) 

-29/-13 

-20/  -7 

-12/  -l 

-4/  4 

3/  10 

10/ 

15 

16/  21 

21/ 

27 

4.0 

{  200.) 

-29/-13 

-21/  -7 

-12/  -2 

-5/  4 

3/  10 

10/ 

15 

16/  21 

22/ 

27 

4.0 

(  600.) 

—30/— 12 

-23/  -7 

-14/  -1 

-5/  3 

1/  10 

9/ 

16 

16/  20 

22/ 

26 

4.0 

(1200.) 

—2  8 /— 1 1 

-20/  -7 

-14/  -1 

-6/  5 

2/  10 

9/ 

15 

15/  21 

21/ 

26 

3.0 

(  50.) 

-23/  -9 

-16/  -3 

-8/  l 

-1/  6 

6/  12 

11/ 

17 

17/  23 

23/ 

28 

3.0 

(  200.) 

-22/-U 

-16/  -4 

-9/  1 

-1/  6 

5/  12 

11/ 

17 

17/  22 

23/ 

28 

3.0 

(  600.) 

-25/  -9 

-17/  -4 

-10/  1 

-2/  6 

4/  12 

10/ 

17 

16/  22 

23/ 

27 

3.0 

(1200.) 

-22/  -8 

-15/  -4 

-10/  2 

-3/  6 

3/  13 

10/ 

18 

16/  23 

22/ 

28 

2.0 

(  50.) 

-17/  -7 

-10/  -l 

-3/  3 

2/  9 

8/  14 

14/ 

19 

20/  24 

25/ 

29 

2.0 

(  200.) 

-17/  -7 

-10/  -2 

-4/  3 

2/  8 

8/  14 

14/ 

19 

19/  24 

25/ 

29 

2.0 

(  600.) 

-18/  -7 

-12/  -2 

-6/  3 

1/  7 

7/  13 

13/ 

19 

19/  24 

25/ 

29 

2.0 

(1200.) 

-17/  -6 

-11/  -1 

-5/  3 

0/  9 

7/  14 

13/ 

19 

19/  24 

24/ 

29 

1.0 

(  50.) 

-10/  -4 

-4/  0 

1/  5 

6/  10 

12/  15 

17/ 

20 

22/  25 

27/ 

30 

1.0 

(  200.) 

-12/  -4 

-6/  0 

0/  5 

6/  10 

11/  15 

17/ 

20 

22/  25 

27/ 

30 

1.0 

(  600.) 

-13/  -5 

-6/  0 

-1/  5 

5/  10 

10/  15 

16/ 

21 

21/  25 

27/ 

30 

1.0 

(1200. ) 

«**/«** 

-6/  0 

-2/  6 

4/  10 

10/  15 

15/ 

20 

21/  25 

26/ 

30 

-36- 


terminated.  In  addition  to  looking  at  the  temperature  of  the  first  echoes, 
the  height  of  the  echoes  was  investigated.  Table  8  shows  the  same  results 
as  Table  7,  but  now  presents  the  echo  heights  in  terms  of  height  above  cloud 
base  (km).  In  this  presentation,  the  effect  of  cloud  base  temperature  on  echo 
height  is  weakened  considerably  as  compared  to  the  parcel  model  previously 
fe*. ployed.  The  water  contents  within  the  cloud  are  lower  when  the  cloud  base 
temperatures  are  reduced  and  the  echo  takes  longer  to  form,  but  the  actual 
position  of  the  echo  that  eventually  does  form  is  not  dramatically  different. 

For  example,  for  the  case  of  a  8  m/sec  updraft  (600  droplets  cm  )  the  parcel 
model  predicted  that  you  might  produce  a  10  dBZ  first  echo  by  3.2  km  above 
cloud  base  for  a  30°C  cloud  base  temperature  and  by  4.0  km  above  cloud  base 
for  a  10°C  base  temperature.  With  the  trajectory  model,  the  base  of  the  echo 
for  the  30°C  cloud  base  temperature  would  be  at  2.1  km  and  top  of  the  echo  at 
3.9  km.  For  a  10°C  cloud  base  temperature,  the  echo  base  would  be  at  2.6  km 
and  the  top  at  4.3  km.  In  this  particular  example  the  sensitivity  to  variations 
in  cloud  base  temperature  is  essentially  cut  in  half. 

Table  9  shows  the  effect  of  reducing  the  liquid  water  contents  to  half 
their  adiabatic  values.  As  would  be  expected,  the  reduction  in  liquid  water 
content  results  in  higher  first  echoes ,  but  the  effect  is  not  as  dramatic  as 
was  predicted  with  the  parcel  model.  Table  10  shews  the  effect  of  further 
reducing  the  water  contents.  In  this  case,  they  have  been  reduced  to  a  function 
of  height  above  cloud  base  specified  by  (1).  As  discussed  above  (see  Background 
Section),  this  expression  gives  an  estimate  of  the  mean  water  content  expected 
in  a  convective  cloud.  As  expected,  echoes  are  higher  still  and  the  conditions 
for  which  the  model  fails  to  predict  echo  formation  are  enlarged. 


TABLE  8 


AX)  I  BAT  >1.0 

UPDRAFT  VELOCITY ,  M/S 
(I  CONCENTRATION,  CM-3)  -5  0 


16.0 

(  50.) 

*** /*** 

***/*** 

16.0 

(  200.) 

**•/**• 

***/««* 

16.0 

(  600.) 

*•*/•** 

•**  / 1  *  « 

16.0 

(1200.) 

***/•*• 

*«*/**« 

14.0 

(  50.) 

***/••* 

5. 5/5. 3 

14.0 

(  200.) 

ft**/*** 

5. 5/5. 4 

14.0 

(  600.) 

ft**/*** 

*»« /»•• 

14.0 

(1200.) 

ft**/*** 

»•*/*•* 

12.0 

(  50.) 

ft**/*** 

5. 4/4. 4 

12.0 

(  200.) 

***/*** 

5. 6/4. 6 

12.0 

(  600.) 

ft**/*** 

5. 4/4. 7 

12.0 

(1200.) 

••*/«** 

5. 4/5.0 

10.0 

(  50.) 

4. 4/3. 6 

5. 5/3. 5 

10.0 

(  200.) 

4. 4/3. 9 

5. 5/3.5 

10.0 

(  600.) 

4. 2/4.0 

5. 4/3. 9 

10.0 

(1200.) 

***/•*« 

5. 4/4.1 

8.0 

(  50.) 

4. 6/3.1 

4. 9/3.0 

8.0 

(  2 CO.) 

4. 4/3. 2 

5. 1/2. 7 

8.0 

(  600.) 

4. 4/3. 3 

5. 1/2. 9 

8.0 

(1200.) 

4. 4/3. 5 

5. 1/3. 2 

6.0 

(  50.) 

4. 3/2. 5 

4. 1/2.2 

6.0 

(  200.) 

4. 4/2. 4 

4. 0/2. 3 

6.0 

(  600.) 

4. 3/2. 2 

4. 3/2. 2 

6.0 

(1200.) 

4. 2/2. 3 

4. 3/2. 2 

4.0 

(  50.) 

3. 3/1.2 

3. 1/1.2 

4.0 

(  200.) 

3. 3/1. 3 

3. 3/1.2 

4.0 

(  600.) 

3. 4/1.1 

3. 4/1.1 

4.0 

(1200.) 

3. 2/1.0 

3. 1/1.2 

3.0 

(  50.) 

2.6/0. 6 

2. 5/0. 6 

3.0 

(  200.) 

2.5/0. 9 

2. 5/0. 7 

3.0 

(  600.) 

2. 8/0. 6 

2. 7/0. 7 

3.0 

(1200.) 

2. 5/0. 4 

2. 5/0. 7 

2.0 

(  50.) 

1.7/0. 3 

1. 7/0.2 

2.0 

(  200.) 

1.8/0. 3 

1.7/0. 3 

2.0 

(  600.) 

1.9/0. 3 

1.9/0. 4 

2.0 

(1200.) 

1.8/0. 2 

1.9/0. 2 

1.0 

(  50.) 

0.8/-. 2 

0.8/-.1 

1.0 

(  200.) 

1.0/-. 1 

1.0/-. 1 

1.0 

(  600.) 

1. 2/0.0 

1. 1/0.0 

1.0 

(1200.) 

»«•/•** 

1. 0/0.0 

10  DBZ  FIRST  ECHO  HEIGHT 

CLOUS  BASE  TEMPERATURE 


5 

10 

15 

20 

25 

7. 0/5. 4 

7. 3/4. 8 

7. 0/4. 6 

6. 5/4. 4 

6. 1/4. 2 

6. 9/5. 4 

6. 8/4. 8 

6. 5/4. 6 

6. 1/4. 2 

6. 2/4.0 

6. 9/5. 6 

7. 2/5.0 

6. 5/4. 6 

6. 1/4. 4 

6. 2/4. 2 

6. 8/5. 8 

7. 1/5. 2 

7. 0/4. 8 

6. 6/4. 6 

6. 2/4.4 

6.6/4. 7 

6. 3/4. 3 

6. 0/4.1 

5. 7/3. 9 

5. 7/3.6 

6. 7/4. 7 

6. 4/4. 2 

6. 1/4.0 

5. 7/3. 8 

S.3/3.6 

6. 5/5.0 

6. 3/4. 4 

6. 0/4. 2 

5. 7/4.0 

5. 8/3. 8 

6. 6/5. 2 

6. 7/4.6 

6. 0/4. 4 

6. 1/4.0 

5. 7/3. 8 

6. 1/4.0 

5. 8/3. 8 

5. 5/3. 5 

5. 2/3. 4 

5. 2/3. 2 

5. 9/4.1 

5. 8/3. 7 

5. 6/3. 4 

5. 2/3. 3 

4. 9/3. 2 

6. 1/4. 2 

5. 7/3. 8 

5. 5/3. 6 

5. 2/3. 4 

4. 8/3.4 

6. 2/4. 5 

6. 0/4.0 

5. 4/3. 8 

5. 1/3.6 

5. 2/3. 4 

5. 5/3. 2 

5. 0/3.1 

4. 8/2. 9 

4. 5/2. 8 

4. 2/2. 7 

5.5/3. 3 

5. 1/3.1 

4. 8/2. 9 

4. 6/2. 8 

4. 3/2. 7 

5. 4/3.5 

4. 9/3. 2 

4. 8/3.0 

4. 5/2. 9 

4. 3/2.8 

5. 5/3. 7 

5. 2/3.4 

5. 1/3. 2 

4. 8/3.0 

4. 5/2. 9 

4. 7/2. 6 

4. 4/2. 4 

4. 2/2. 3 

4. 0/2. 2 

3. 8/2. 2 

4. 7/2. 5 

4. 4/2. 4 

4. 0/2. 4 

3. 8/2. 3 

3. 8/2.1 

4. 6/2. 9 

4. 3/2.6 

4. 2/2. 4 

4. 0/2. 3 

3.B/2.2 

5. 0/2. 9 

4. 4/2. a 

4. 4/2.5 

4. 2/2. 4 

4.0/2. 3 

3. 7/2.1 

3. 4/2.0 

3. 3/1.8 

3. 3/1. 7 

3. 1/1. 7 

3. 7/2.1 

3. 5/2.0 

3. 4/1. 7 

3.3/1. 7 

3.1/1. 6 

4. 0/1. 9 

3. 6/1. 9 

3.5/1. 7 

3.2/1. 8 

3. 3/1.6 

4. 0/2.1 

3. 9/1. 9 

3. 7/1. 8 

3. 3/1.9 

3. 2/1.8 

3. 0/1.1 

2. 7/1.3 

2. 6/1.0 

2.4/1. 2 

2. 3/1. 2 

2. 9/1. 4 

2. 8/1.1 

2.6/1. 1 

2. 4/1. 2 

2. 4/1. 2 

3. 3/1.1 

2. 9/1. 4 

2. 7/1.2 

2. 7/1.1 

2. 5/1. 2 

3.2/1. 1 

3. 1/1.0 

2. 8/1.1 

2. 7/1.2 

2. 7/1.1 

2.4/0. 7 

2. 2/0. 8 

2. 1/0. 7 

2. 1/0. 7 

2. 1/0.6 

2.5/0. 7 

2.3/0. 7 

2.2/0. 7 

2. 1/0. 7 

2. 1/0. 7 

2.6/0. 7 

2. 4/0. 9 

2. 4/0. 7 

2. 3/0. 7 

2. 3/0. 7 

2. 6/0. 6 

2. 5/0. 7 

2. 6/0. 5 

2. 5/0. 5 

2. 4/0. 6 

1.5/0. 3 

1.5/0. 3 

1. 5/0.2 

1.5/0. 2 

1.4/0. 3 

1.7/0. 3 

1.6/0. 4 

1. 6/0.3 

1.6/0. 2 

1.5/0. 3 

1. 9/0.4 

1. 7/0.5 

1. 7/0,4 

1. 7/0.3 

1. 7/0.3 

1.9/0. 3 

1.9/0. 3 

1.8/0. 3 

1. 7/0.4 

1. 7/0.4 

0.8/-.1 

0. 6/0.0 

0. 8/0.0 

0.8/-. I 

0. 8/0.0 

0. 9/0.0 

0. 9/0.0 

0. 8/0.0 

0. 9/0.0 

0. 8/0.0 

1. 1/0.0 

1. 0/0.1 

1. 1/0.1 

1.1/-. 1 

1. 0/0.0 

1.3/-. 2 

1. 3/0.0 

1. 2/0.0 

1. 2/0.0 

l. 2/0.0 

(TOP/BASE) 


iP 

ti. 1/4.0 
S. 7/4.0 
5. 7/4.0 
6. 2/4. 2 

S.3/J.7 
5. 4/3. 6 
5. 4/3.6 
5. 3/3.8 

4. 9/3. 2 
4. 9/3.0, 
4, 9/3. 2 
4. 9/3. 4 

4. 3/2. 6 
4. 4/2. 6 
4. 3/2. 7 
4. 3/2. 8 

3. 9/2.0 
3. 6/2.1 
3. 9/2.1 
3.8/2. 3 

3. 1/1.5 
2.9/1. 7 
3. 1/1.6 
3. 3/1.6 

2.4/0. 9 
2. 3/1. 2 
2. 4/1. 2 
2. 6/1.1 

2. 0/0. 6 
2. 1/0. 7 
2. 1/0. 8 
2. 4/0. 5 

1. 4/0.3 
1.5/0. 3 
1.6/0. 4 
1.8/0. 2 

0. 8/0.0 
0. 8/0.0 
0. 9/0.0 
l. 2/0.0 


TABLE  9 


ADI BAT  -  0.5 


10  DBZ  FIRST  ECHO  KEICHT  (TOF/BASE) 


CLOUD  BASE  TEMPERATURE 


UPDRAFT  VELOCITY,  M/S 


(1  CONCENTRATION,  CM-3) 

-5 

0 

5 

10 

15 

20 

25 

30 

16.0 

(  SO.) 

***/*** 

**•/**« 

7. 7/7.4 

9. 1/6. 8 

8. 9/6. 3 

8. 7/6.0 

8. 2/5. 8 

16.0 

(  200.) 

***/*** 

•••/*•• 

***/«*• 

7. 8/7. 4 

8. 9/6.7 

8.5/6. 3 

8. 2/6.0 

7. 8/5. 9 

16.0 

(  600.) 

**• /*•* 

***/••* 

***/*** 

***/*•« 

9. 0/7.0 

8. 9/6. 6 

8. 6/6. 2 

8. 3/6.0 

16.0 

(1200.) 

**•/**• 

•**/••• 

*•*/**« 

*«*/*«* 

9. 0/7. 4 

9. 2/6. 8 

8. 5/6.4 

8. 1/6. 2 

14.0 

(  50.) 

***/*** 

*•*/*•• 

***/*•• 

7. 8/6. 6 

8. 3/6.0 

8. 2/5. 6 

7. 9/5.4 

7. 6/5. 2 

14.0 

(  200.) 

***/•*• 

**•/••• 

***/•*• 

7. 8/6. 5 

8.475.9 

7. 8/5. 7 

7. 5/5. 3 

7. 2/5. 2 

14.0 

(  600.) 

•*«/*** 

***/*** 

«••/»•* 

7.776.8 

8. 6/6. 2 

8. 1/5. 8 

7. 9/5.6 

7.675.4 

14.0 

(1200.) 

***/*** 

***/••* 

»•*/*«* 

7. 7/7. 3 

8. 7/6. 6 

8.476.2 

7. 7/5.  a 

7.4/5. 6 

12.0 

(  50.) 

•**/**• 

•••/•*• 

6. 4/5. 6 

7. 7/5. 5 

7. 4/5. 2 

7. 4/4. 8 

7. 0/4. 5 

6. 8/4. 5 

12.0 

(  200.) 

***/*** 

*** /*»* 

6. 4/5. 8 

7.675.5 

7. 5/5.1 

7. 0/4. 9 

6.B/4.7 

6.4/4. 5 

12.0 

(  600.) 

***/*«« 

*•• /«•* 

6. 2/6.1 

7. 6/5. 8 

7. 7/5. 3 

?. 2/5.0 

7. 0/4. 8 

6. 8/4. 6 

12.0 

(1200.) 

***/*«• 

*•*/•*• 

*••/**• 

7. 6/6.1 

7. 8/5. 7 

7.4/5. 3 

7 .4/5.1 

7.0/5. 0 

10.0 

(  50.) 

***/*** 

4.4/3. 3 

6. 4/5.0 

7. 0/4. 6 

6.774.2 

6. 3/4.0 

6. 1/3. 9 

5. 8/3.8 

10.0 

(  2 00.) 

*»*/**• 

••*/*•• 

6. 4/4. 8 

7. 0/4. 6 

6. 7/4.2 

6. 3/4.1 

6. 1/3. 9 

5. 9/3. 7 

10.0 

(  600.) 

•** /*•• 

6. 4/5. 2 

7. 2/4. 6 

6. 9/4. 4 

6. 5/4. 2 

6. 0/4.1 

5. 7/4.0 

10.0 

(1200.) 

***/*** 

***/*** 

6. 4/5. 5 

7. 1/5. 2 

6. 9/4. 7 

6. 6/4. 4 

6. 1/4. 4 

6. 0/4. 2 

8.0 

(  50.) 

*«*/•** 

5. 0/3. 7 

6. 1/4.1 

6. 2/3. 4 

5. 7/3.4 

5. 3/3.4 

5. 2/3. 2 

5.0/3. 1 

8.0 

(  200.) 

*** /*** 

5. 0/3. 8 

6. 2/3. 9 

6. 1/3. 6 

5. 7/3.5 

5. 3/3. 4 

5. 2/3.2 

5.0/3. 1 

8.0 

(  600.) 

***/*** 

4. 6/3. 9 

6. 1/4. 2 

6. 0/4.0 

5. 8/3. 7 

5. 5/3. 5 

5. 4/3. 3 

5. 2/3. 2 

8.0 

(1200.) 

***/*** 

4. 2/3. 7 

6. 3/4. 2 

6. 1/4.1 

6. 0/3. 8 

5. 8/3. 5 

5. 4/3. 5 

5. 3/3. 4 

6.0 

(  SO.) 

2. 8/1. 9 

5. 1/3.0 

5. 3/2. 9 

4. 9/2. 8 

4. 5/2. 8 

4. 5/2. 5 

4.272.6 

4. 1/2. 4 

6.0 

(  200.) 

2. 8/2.1 

5. 1/2. 3 

5. 4/2. 9 

5. 1/2. 7 

4. 8/2. 6 

4. 6/2. 6 

4. 2/2. 6 

4. 1/2. 5 

6.0 

(  600.) 

1.8/1. 4 

5. 0/2. 6 

5. 4/2. 9 

5.2/Z.8 

5. 0/2. 6 

4. 5/2. 8 

4. 6/2. 5 

4. 5/2. 3 

6.0 

(1200.) 

•**/**• 

4. 7/2. 6 

5. 3/2.8 

5. 0/3. 2 

5. 0/2. 7 

4. 7/2. 7 

4. 4/2. 7 

4. 4/2. 5 

4.0 

(  50.) 

3. 5/1. 8 

4. 0/1. 6 

3. 9/1. 5 

3. 8/1. 5 

3. 4/1. 8 

3. 4/1. 7 

3. 4/1. 5 

3. 2/1. 6 

4.0 

(  200.) 

3.4/1. 3 

3. 9/1.8 

3. 8/1. 9 

3. 9/1.6 

3. 4/1. 8 

3.5/1. 6 

3. 3/1. 7 

3.3/1. 6 

4.0 

(  600.) 

2. 7/1.1 

3. 9/1. 3 

4. 0/1. 6 

4. 1/1. 4 

3. 9/1. 5 

3. 7/1.6 

3. 7/1.5 

3. 5/1. 6 

4.0 

(1200.) 

2.6/0. » 

3. 5/1. 4 

3. 8/1. 4 

3.5/1. 7 

3. 9/1. 4 

3.8/1. 4 

3. 5/1. 6 

3. 4/1. 6 

3.0 

(  50.) 

3. 2/0. 8 

3.0/1. 1 

3. 0/0. 9 

2. 8/1.1 

2. 8/1.0 

2. 8/0. 9 

2. 7/0.9 

2. 6/1.0 

3.0 

(  200.) 

3. 1/0. 8 

3. 1/1.0 

3. 0/1. 2 

2. 9/1.1 

3. 0/1.0 

2. 8/1.1 

2. 8/1.0 

2. 7/1.1 

3.0 

(  600.) 

2. 6/0. 5 

3. 1/0. 6 

3. 2/1.0 

3. 3/0.9 

3. 1/1.1 

3. 1/1.0 

2. 8/1. 2 

2. 9/1. 2 

3.0 

(1200.) 

2. 3/0. 5 

2.5/0. 7 

2. 7/0. 8 

2. 8/1.0 

3. 0/0. 9 

3. 1/0.8 

3. 0/0. 9 

3. 0/0. 9 

2.0 

(  50.) 

2. 2/0. 2 

2. 2/0. 3 

2. 0/0. 6 

2. 0/0. 4 

1.9/0. 5 

1.9/0. 4 

1.8/0. 5 

1.9/0. 4 

2.0 

(  200.) 

2. 2/0. 3 

2. 2/0. 5 

2. 2/0. 5 

2. 0/0. 6 

2. 1/0. 5 

2. 0/0. 6 

2. 0/0. 5 

2. 0/0. 5 

2.0 

(  600.) 

1.7/0. 2 

2. 0/0. 3 

2. 2/0. 5 

2. 4/0. 5 

2. 3/0.6 

2. 3/0. 5 

2.1/0. 7 

2. 2/0. 6 

2.0 

(1200.) 

***/*** 

2. 0/0. 2 

2. 0/0. 5 

2. 2/0. 3 

2. 1/0. 5 

2. 0/0. 6 

2. 1/0. 5 

2. 1/0. 5 

1.0 

(  50.) 

*••/••* 

**• /•*• 

0.9/-. 2 

0. 9/0.0 

0. 9/0.0 

0. 9/0.0 

0.9/-.1 

0.9/-.1 

1.0 

(  200.) 

*** /**• 

1. 2/0.0 

1. 2/0.0 

1. 2/0.0 

1. 1/0.1 

1.1/0. 2 

1. 1/0.0 

1.0 

(  600.) 

•*•/••* 

*•«/•*• 

•*•/*•* 

*** /*•• 

1.4/0. 2 

1. 4/0.0 

1. 3/0.1 

1.2/0. 2 

1.0 

(1200.) 

*** /*** 

/«•* 

*•*/*•• 

«*•/*** 

**•/**• 

ft**/*** 

1. 3/0.1 

1. 3/0.0 
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UPDRAFT  VELOCITY,  M/S 


(t  CONCENTRATION,  04-3) 

-5 

0 

s 

10 

15 

20 

25 

30 

16.0 

(  50.) 

***/*** 

***/*•» 

ft**/*** 

ft**/*** 

•*•/**• 

***/9.7 

***/9.2 

**»/8.7 

16.0 

(  200.) 

•**/**• 

*•*/**• 

ft**/*** 

*** /*•* 

«**/**• 

«*V9.8 

**»/9.3 

•••/8.7 

16.0 

(  600.) 

••*/*** 

*••/*** 

ft**/*** 

ft**/*** 

***/*•• 

•••/10. 

***/9. 7 

••*/ 9.2 

16.0 

(1200.) 

*** /**• 

*••/*** 

ft**/*** 

ft**/*** 

***/«*« 

>••/••• 

»•»/«•• 

•••/9.8 

14.0 

(  50.) 

**•/*«* 

**•/•** 

***/*** 

ft**/*** 

8. 6/8. 3 

***/8.6 

***/7.9 

***/7.9 

14.0 

(  200.) 

**•/*•* 

**« /*** 

ft**/*** 

*** /*** 

**•/**« 

***/8.4 

***/8.2 

***/7.7 

14.0 

(  600.) 

*•*/*•* 

*••/*** 

ft**/*** 

ft**/*** 

*««/#«* 

***/9.0 

***/8.5 

***/8.2 

14.0 

(1200.) 

***/*** 

**•/•** 

ft**/*** 

ft**/*** 

«**/**• 

•••/ 9.6 

***/9.1 

••*/8.6 

12.0 

(  50.) 

***/*** 

***/•«* 

ft**/*** 

ft**/*** 

8. 8/7.1 

***/7.1 

***/7.0 

**»/6.7 

12.0 

(  200.) 

***/*•• 

*** /*** 

*•*/*•* 

ft**/*** 

8. 8/7. 2 

***/7.2 

•**/6.8 

***/6.8 

12.0 

(  600.) 

**•/*•* 

***/*** 

ft**/*** 

•«*/*** 

8. 6/8.0 

***/7.5 

***/7.1 

•**/7.0 

12.0 

(1200.) 

•*• /*** 

*•*/**• 

***/*** 

ft**/*** 

•*•/•** 

•**/8.1 

***/7.6 

*»*/7.4 

10.0 

(  50.) 

•*«/*** 

**•/*•* 

ft**/*** 

6. 8/5. 6 

8.8/5. 7 

9.0/S. 7 

8. 6/5. 5 

8.2/5.* 

10.0 

(  200.) 

•**/*** 

•**/•** 

ft**/*** 

6. 8/6.0 

8.8/5. 9 

9.2/5. 6 

8. 4/3.7 

8, 2/5.5 

10.0 

(  600.) 

***/•*• 

*••/•** 

ft**/*** 

«*»/«.. 

8. 4/6. 5 

8. 9/6. 2 

9. 1/5. 8 

8. 5/5. 8 

10.0 

(1200.) 

*••/*** 

*•*/•** 

ft**/*** 

8. 6/6. 8 

9. 1/6. 5 

9.0/6. 2 

8. 6/6.1 

8.0 

(  SO.) 

***/*** 

***/*** 

4.4/3. 7 

7. 1/4. 7 

8.0/4. 6 

7.4/4. 8 

6. 9/4. 8 

7.0/4. 2 

8.0 

(  200.) 

*••/•*• 

***/•*• 

4. 0/3.6 

7. 0/4. 6 

7. 9/4. 8 

7. 4/4. 7 

7. 1/4. 6 

6. 9/4. 5 

8.0 

(  600.) 

***/*** 

•**/••* 

*•*/*•* 

6. 8/4. 8 

7. 6/5.1 

7. 9/4. 7 

7. 5/4. 8 

7.0/4. 7 

8.0 

(1200.) 

*•*/*•* 

•*•/*** 

***/*** 

6. 6/5. 2 

7. 9/5.0 

7. 8/4. 9 

7. 7/4. 8 

7.4/4. 7 

6.0 

(  50.) 

**•/**• 

I. 0/0. 4 

5.4/2. 9 

6.4/3. 3 

6.0/3.  7 

5. 8/3. 6 

5. 9/3. 3 

5. 4/3.5 

6.0 

(  200.) 

*••/**« 

0.2/-. I 

5. 0/2. 6 

6. 4/3.1 

6. 5/3. 2 

6. 2/3.4 

6. 0/3. 3 

5. 8/3.1 

6.0 

(  600.) 

••*/**• 

*«•/•** 

4. 2/2. 6 

5.B/3.2 

6. 2/3. 2 

6. 4/3.1 

6. 3/3.3 

6.0/3. 3 

6.0 

(1200.) 

••*/*•• 

*•*/*•• 

4. 0/2. 5 

5. 6/3. 2 

5. 9/3. 3 

6. 1/3.1 

6. 1/3. 2 

6. 0/3. 2 

4.0 

(  50.) 

••*/«*• 

3. 3/1.1 

4. 1/1.9 

4. 4/1. 9 

4. 3/1. 8 

4. 2/1. 8 

4.1/1. 7 

4. 1/1. 6 

4.0 

(  200.) 

•**/••* 

2.6/0. 7 

4. 4/1. 3 

4. 4/1. 9 

4. 5/1. 9 

4. 2/2.0 

4. 2/1.8 

4. 1/2.0 

4.0 

(  600.) 

••*/*•* 

1.8/0. 7 

3. 5/1.0 

4. 0/1. 3 

4. 1/1. 7 

4. 5/1. 5 

4. 4/1. 6 

4. 5/1. 6 

4.0 

(1200.) 

•**/*•• 

*••/*•* 

3. 3/1.2 

3. 6/1. 4 

3. 8/1. 6 

3. 7/1.8 

4. 0/1. 7 

4. 0/1. 6 

3.0 

(  50.) 

*•*/*•* 

2.7/0. 5 

3. 1/0.8 

3. 3/0. 8 

3. 3/1.0 

3. 3/0. 9 

3. 1/1.1 

3. 2/0.9 

3.0 

(  200.) 

**•/*•* 

2. 4/C. 4 

3.2/0.? 

3.4/1. 1 

3. 2/1. 2 

3. 3/1.1 

3. 2/1.1 

3. 1/1. 2 

3.0 

(  600.) 

*•*/•** 

•••/••* 

2. 3/0. 3 

2. 7/0.6 

3.1/0. 7 

3.3/0. 7 

3. 2/0. 8 

3. 2/1.1 

3.0 

(1200.) 

•*•/*•* 

•••/**• 

2. 2/0. 6 

2. 5/0. 6 

3. 0/0. 5 

3. 1/0. 6 

3.0/0. 7 

3.0/0. 7 

2.0 

(  50.) 

**•/**• 

l. 4/0.1 

2. 0/0.0 

2. 1/0.1 

2. 2/0.3 

2.0/0. 3 

2.0/0. 3 

2.0/0. 2 

2.0 

(  200.) 

*•*/••• 

•*•/•*« 

1. 7/0.3 

2. 0/0. 3 

2. 0/0. 5 

2. 2/0. 3 

2. 1/0. 4 

2. 1/0. 5 

2.0 

(  600.) 

•ft*/*** 

*•*/**• 

•**/••* 

1. 7/0.1 

2. 0/0.0 

1.9/0. 3 

2. 1/0.3 

2. 2/0.1 

2.0 

(1200.) 

•**/••* 

**•/*•• 

**«/««. 

*«*/«** 

1.8/0. 2 

2.0/0. 3 

1.9/0. 2 

1.9/0. 2 

1.0 

(  50.) 

•••/•*• 

•**/•** 

***/*** 

**•/•*• 

**«/**• 

0.6/-.1 

0.7/-. 3 

1.0 

(  200.) 

*•• /**• 

••*/*** 

***/*•• 

••*/*** 

•*•/**• 

»**/•** 

1.0/-. 2 

1.0 

(  600.) 

•**/*•• 

•••/*** 

*••/*** 

***/•*• 

*** /••* 

**• /*** 

*•* / *  *  * 

1.0 

(1200.) 

***/*•• 

**•/**• 

•«*/ ... 

*«*/*•* 

**•/*•* 

tM/m 

*••/«•* 

»«*/*** 

One  of  the  surprising  results  that  was  discovered  in  this  analysis  was 
the  potential  for  growth  of  large  liquid  drops  in  relatively  cold-based 
convective  clouds.  Such  clouds  are  typical  of  the  High  Plains  (Colorado, 

Montana,  etc.  ).  Although  precipitation  formation  in  these  clouds  is  primarily 
an  ice  or  mixed-phase  process,  large  liquid  drops  have  been  observed  in  these 
•*louds.  Since  it  was  thought  to  be  impossible  for  such  drops  to  grow  directly 
by  coalescence,  these  observations  were  interpreted  as  evidence  for  recycling  of 
ice  or  graupel  that  fall  below  the  0°  isotherm,  melt,  and  are  then  carried  back 
into  the  updraft.  The  current  model  results,  on  the  other  hand,  suggest  that 
direct  growth  of  liquid  drops  is  not  only  possible,  but  likely. 

Figure  17  illustrates  some  of  the  large  drop  trajectories  in  the  case 

of  a  5  m/sec  updraft  and  a  5°C  cloud  base  temperature.  Figure  18  shows,  for 

this  run,  the  evolution  of  the  expected  size  distribution  of  liquid  drops 

that  would  be  found  at  the  -5°C  and  -10°C  levels  of  the  cloud.  The  small  drops 

are  the  first  to  arrive  at  a  given  level,  with  the  larger  drops  lagging  behind. 

Eventually  the  concentration  of  large  drops  builds  up,  as  can  also  be  seen  in 

-3 

Figure  19  which  shows  the  concentrations  (m  }  of  half  millimeter  (left  hand 
figure)  and  millimeter  (right  hand  figure)  diameter  drops.  Even  at  higher 
updraft  velocities  it  is  possible  to  grow  large  liquid  drops  during  the  primary 
ascent.  Figure  20,  for  example,  shows  the  possible  concentration  of  millimeter 
drops  (top  figures)  and  half  millimeter  drops  (bottom  figures)  -or  updrafts  of 
10  m/sec  (left  hand  figures),  15  m/sec  (middle  figures),  and  20  m/sec  (right  hand 
figures).  In  each  case,  the  initial  concentration  of  large  droplets  that  grow 
to  raindrop  dimensions  was  based  on  aerosol  observations  that  were  made  over  the 
High  Plains  during  the  HIPLEX  project. 
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Figure  17.  Sample 
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Figure  18.  Evolution  of  the  drop-size  distribution  as  a  function  of  time 
and  height  above  cloud  base. 
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ure  19.  Concentration  (m  )  of  half  millimeter  (left 'figure)  and 
millimeter  (right  figure)  diameter  drops. 
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:  20.  Steady-state  concentration  of  millimeter  (top  figures)  and  half 

millimeter  (bottom  figures)  drops  in  10  (left  figures),  15  (center) 
and  20  m/sec  (right)  updrafts. 
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ICE-MULTI PLICATION 


In  some  clouds,  many  more  ice  crystals  are  observed  at  temperatures 
of  -10°C,  or  warmer,  than  can  be  accounted  for  by  measured  concentrations 
of  ice  nuclei.  This  has  long  been  interpreted  as  evidence  for  some  sort  of 
ice  multiplication  mechanism.  In  a  series  of  articles,  Hallett  and  Mossop 
have  reported  on  laboratory  experiments  supporting  a  possible  "splintering" 
mechnaism  for  the  production  of  secondary  ice  crystals  (Hallett  and  Mossop, 
1974;  Mossop  and  Hallett,  1974;  Mossop,  1978a;  Mossop  and  Wishart,  1978). 

In  laboratory  studies  of  ice-multiplication  Mossop  (1978a)  suggested 

that  the  rate  of  production  of  secondary  ice  crystals  is  proportional  to 
0  93 

N^N's  ’  ,  where  is  the  number  of  large  drops  (D  >  24  ym)  swept  up 

each  second  by  a  falling  graupel  particle  at  -5°C  and  N’s  is  the  corresponding 
sweep-out  rate  of  small  (D  <  13  ym)  drops.  Recently  Mossop  (personal 
communication)  has  suggested  that  the  exponent  Q.93  may  be  unrealistically 
precise  and  that  a  simplier  proportionality  of  N  is  probably  more 
appropriate . 

In  many  cases,  the  cloud  droplet  distribution  can  be  descirbed 
accurately  by  a  gamma  distribution.  Following  Berry  and  Reinhardt  (1974), 
for  example,  the  cloud  droplet  distribution  can  be  expressed  as 

f<x>  =  j-  G(v)  S'v  e'(1*v5S  (5) 

where. 


G(v) 


(1+v) 


hv 


r (i+v) 


and  where  f<x>  is  the  number  density  of  drops  of  mass  x,  N  is  the  total 
number  concentration  of  cloud  droplets,  L  is  the  liquid  water  content, 
and  v  is  a  shape  parameter  specifying  the  breadth  of  the  distribution. 

The  relative  dispersion  (o  /r)  of  droplet  distribution  is  often  observed 
to  be  near  0.2.  This  corresponds  rather  closely  to  a  value  of  v  =  2. 

In  some  observations,  however,  the  relative  dispersion  appears  to  increase 
with  height  above  cloud  base  (e.g.,  Warner,  1969).  In  this  case  the 
appropriate  value  of  v  could  vary  from  v  =  2  down  to  v  =  0  (relative 
dispersion  =  0.4)  . 

Whatever  the  value  of  v,  of  course,  it  is  always  possible  to  integrate 
(5)  to  give  the  number  of  drops  larger  or  smaller  than  any  specified  size 
limit.  This  means  that,  with  an  appropriate  choice  for  v,  the  total  droplet 
distribution  can  be  expressed  in  terms  two  physically  meaningful  parameters: 

N  and  L.  The  relative  number  of  large  and  small  drops  swept  out  by  a 
falling  graupel  particle,  however,  will  not  necessarily  be  directly 
proportional  to  the  relative  number  concentration  of  large  and  small  drops 
since  the  collection  efficiency  for  capture  by  the  graupel  will  be  different 
for  the  large  and  small  size-ranges  that  are  of  interest.  Examination  of 
the  collection  efficiencies  suggest  that  the  collection  efficiencies  for 
small  drops  ( D  <  1 3  pm)  is  typically  about  80%  of  the  corresponding  collection 
efficiency  for  large  drops  (D  >  24  pm)  over  a  wide  range  of  possible  gTaupel 
diameters . 

Taking  this  relative  collection  efficiency  into  account,  figures  21, 

13,  s.vi  25  illustrate  the  relative  effectiveness  of  the  Hal lett-Mossop  ice- 
multiplication  mechanism  as  a  function  of  the  cloud  base  temperature  and 
droplet  concentration.  In  all  cases  the  contours  of  relative  efficiency 
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Relative  efficiency  of  the  Hallett-Mossop 
i ce -mul t ip  1 icat ion  mechanism  for 
adiabatic  water  contents. 


CLOUD  BASE  DROPLET  CONCENTRATION  (cm- 


CLOUD  BASE  DROPLET  CONCENTRATION  (cm~ 


50%  of  ADIABATIC 


CLOUD  BASE  TEMPERATURE  (°C) 


Relative  efficiency  of  the  Hal  1 ett-Mossop 
ice-multiplication  mechanism  for  water 
contents  that  are  50°o  of  adiabatic. 
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are  labeled  in  relative  units  only.  Figure  21  is  based  on  adiabatic  water 
contents,  while  Figure  23  is  appropriate  for  water  contents  that  are  half 
adiabatic  and  Figure  25  corresponds  to  the  mean  observed  water  contents 
shown  in  Figure  2.  In  each  of  these  cases,  the  shape  parameter  v  was 
assumed  to  be  constant  (=2)  .  Figures  22,  24,  and  26  repeat  the  same  sort 
of  estimates  of  ice-multiplication  effectiveness,  but  in  this  case  invoke 
a  variable  v  which  ranges  from  2  (at  cloud  base)  to  1  (at  750  m  above 
cloud  base)  and  ultimately  to  0  (at  2  km  altitude  and  above).  These  plots 
agree  rather  well  with  the  empirical  "boundary”  .suggested  by  Mossop  (l&78b) 
separating  conditions  conducive  for  ice  multiplication  from  those  not 
conducive,  but  goes  well  beyond  the  boundary  concept. 


SNOWFLAKE  AGGREGATION 


In  most  detailed  microphysical  modeling  efforts,  snowflake  aggregation  is 

simply  ignored.  In  recent  years,  however,  aggregation  has  increasingly  been 

spotlighted  as  an  important  physical  process  that  needs  to  be  better  understood. 

While  in  principle,  aggregation  of  snowflakes  into  larger  and  larger  units  seems 

to  be  a  rather  straightforward  problem,  the  wide  variety  of  particle  shapes, 

bulk  densities,  and  fall  behavior  as  well  as  the  potential  interactions  between 

nearby  particles  have  presented  formidable  obstacles  to  anything  more  than  a 

qualitative  understanding  of  the  problem.  As  recently  as  1971,  Mason  said, 

"The  growth  of  these  aggregates  is  governed  by  the  collision  and 
aggregation  efficiencies  of  the  crystals  and  by  their  relative 
motions  which,  because  of  the  aerodynamic  problems  posed  by  the 
complexity  and  variability  of  the  ice-crystal  geometry,  are  not 
amenable  to  quantitative  computation." 

Rassarelli  (1978a),  has  developed  an  analytical  approach  to  modeling  the 
evolution  of  the  snow-flake  size  distribution  by  vapor  deposition  and  aggregation 
of  flakes  chat  bypasses  many  of  the  difficulties  mentioned  by  Mason.  This 
approach  involves  assuming  an  exponential  form  for  the  snowflake  size  distri¬ 
bution, 


_3 

where  n(D)dD  is  the  number  of  snowflakes  (cm  )  in  the  diameter  range  from 
(D-50D)  to  (D+ldD),  and  then  solving  the  moment  conservation  equations  for  the 
total  mass  and  radar  reflectivity  in  terms  of  the  distributional  parameters 
Ne  and  .  The  end  result  is  to  produce  expressions  for  these  parameters  as 
functions  of  time.  Although  this  approach  ignores  many  of  the  complexities  that 
can  complicate  the  analysis  (e.g.  Rassarelli  and  Srivastava, 


1879),  Passarelli 
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has  had  considerable  success  in  applying  this  analytical  approach  to  observational 
data  (e.g.  Passarelli,  1978b;  Lo  and  Passarelli,  1982) 

If  the  snowflake  sice  distribution  is  truely  exponential,  then  the 
distributional  parameters  N0  and  X  can  be  related  to  other  more  physically 
meaningful  properties  of  the  distribution  such  as  the  total  number  (n^)  or 
mass  (m^)  of  snowflakes  per  unit  volume. 


it  pL  Nc 


where  is  the  snowflake  bulk  density.  Alternately,  the  number  or  mass  of  snow¬ 
flakes  larger  than  a  given  diameter,  D0 ,  can  be  expressed  as  : 


n ( D > D_  ) 


m ( D  >  C0  ) 


)r,  D02  X2  D  3  X  3 

it.  e  1  +  D  0  X  +  -  +  - 

2  6 


In  r.anv  cases,  these  last  two  parameters  are  cf  mere  interest  chan  (~)  and  (■). 
In  terms  of  precipitation  formation,  the  parameter  of  interest  is  the  rate  cf 
increase  of  the  mass  cf  large  flakes  which  can  be  expressed  as 


dm  (  D  >  Dr  ) 

dt  ~ 


-h  D“  X3  e_XD°  — 

T  0  dt 


’itical  term  in  <  '  1  ) is  dX/dt,  which  can  be  obtained  from  Passarelli  (1978a). 


For  the  case  of  aggregation  alone, 


where  b  is  the  exponent  in  a  particle  fallspeed  relation  of  the  form 


v 

t 


The  constant  in  (12)  is  defined  as 


a  D 


(13) 


_  (1-b)  a  E  Kb) 

2  "  2880  p 

where  a  and  b  are  defined  by  the  fallspeed  relation  (13 ),  E  is 
efficiency,  and  is  the  snowflake  bulk  density.  Both  E  and 
be  independent  of  particle  size.  The  term  1(b)  is  defined  by 
relation: 


(iu) 

a  mean  collection 
p .  are  assumed  to 

l 

a  rather  intimidating 


(b) 


(p)  2 


1-p 


a 

y  c 

F ( 1  ,p  ;  8-i ; 5 ) 

F( 1  ,p  ;  4+b+i ; 5 ) 

1 .  i 

—  X 

L  7-i 

3+b  +  i 

(15) 


where  p=10+b  ,  C_ =C  =1  ,  C  =2  ,  T  represents  the  gamma  function,  and  F  represents 
Gauss'  hypergeometric  function.  For  these  arguments,  the  series  solution  for  F 
can  be  expressed  as 


F  ( 1 


;  c  ;  2  ) 


r_(c  )  y  r  (b+n  )  ,  ^  n 
r (b )  r ( c  +n  ) 


( i  -: ) 


While  (15)  and  (15)  can  be  readily  evaluated  on  modern  computers,  the  constant 
1(b)  itself  is  rather  well  behaved  ever  the  range  of  interest  (see  Figure  27) 
and  can  also  be  estimated  using  a  simple  polynomial  expression: 


1(b)  =  a1b+a2b2+a„b3 


where  a.  =  1.57  x  10' 


=  1.86  x  10' 


and  a. 


2.87 


(”) 
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r ig ure  53  shows  an  example  of  a  series  of  observations  of  snowflake  size 
distributions  that  were  obtained  in  a  advecting  spiral  descent  (Lo  and  Passarelli 
1982),  plotted  as  functions  of  the  individual  slope,  X,  and  intercept,  N0 , 
parameters.  The  general  trend  of  the  observations  during  the  descent  is  from 
t;e  upper  right  (highest  observations)  to  lower  left  (the  lowest  heights 
observed).  Figure  29  shows,  as  functions  of  the  same  distributional  parameters, 
tne  rate  of  increase  of  large  snowflake  mass  (g  cm  sec  )  for  the  case  E  =  1.4, 
a  =  150  (cgs),  b  =  0.31,  Dc  =  0.3  cm,  and  p^  =  0.05  g  cm  .  For  the  same  case. 
Figure  30  shows  the  proportional  rate  of  increase  of  snowflake  mass  in  large 
flakes 

1  dm( D>Dn ) 
m(D>D0 )  dt 

defined  by(10)ar.d  (11),  which  is  often  a  mere  useful  way  of  expressing  the 
transfer  of  mass  to  large  particles  and  was  extensively  used  in  the  warm  rain 
initiation  aspects  of  this  study. 

,-r.  interesting  way  to  cast  figures  of  this  sort  in  a  new  light  is  to  convert 
them  to  a  new  coordinate  system  that  moves  away  from  the  slope  and  intercept 
car ureters  in  favor  cr  more  physically  meaningful  properties  such  as  number  or 
mass  concentrations  of  crystals.  Figures  31  and  32,  for  example,  illustrate  the 
same  relationships  as  were  shown  in  figures  29  and  30,  but  in  terms  of  number  and 
mass  rather  than  abstract  "parameters."  In  this  case,  the  number  and  mass  were 
define:  using  3  ar.i  10  (-'o'  1.C3  cm),  which  is  a  convenient  threshold  for  data 
sets  obtained  from  airborne  tarticle  probes.  Most  interestingly,  this  conversion 
susses ts  that  it  may  be  possible  to  parameterize  the  growth  of  large  aggregates 
as  a  relatively  simple  function  of  the  total  ice  particle  mass —  the  sane  concept 
u:ed  in  Kessler-style  p ar -mater izat ions  of  warm-cloud  microphysics. 


.v 


,  . 


Alternately , 


Parameter,  cm 


a  illustrating  snowflake  aggregation 
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Abstract 


Riming  growth  rates  for  graupel  and  frozen  raindrops  are  compared 
with  those  of  unfrozen  drops  growing  by  coalescence.  The  results  indi¬ 
cate  that  graupel  can  have  a  considerable  advantage  over  unfrozen 
raindrops  in  regions  where  the  cloud  droplets  are  large.  On  the  other 
hand,  if  the  cloud  droplets  are  small  the  situation  can  reverse  itself 
with  water  drops  being  favored  over  graupel.  In  addition,  the  results 
suggest  that  low  density  graupel  or  aggregates  may  be  particularly 
effective  growth  centers  for  accretion,  and  that  particles  entering  a 
wet  growth  regime  may  experience  a  sudden  jump  in  growth  associated  with 
a  change  in  the  particle's  surface  roughness. 


1.  Introduction 


More  than  20  years  ago,  Braham's  (1954)  observations  of  precipitation 
development  in  Missouri  clouds  led  him  to  conclude  that  snow  pellets 
(graupel)  grew  faster  by  riming  than  unfrozen  drops  grew  through 
coalescence.  This  conclusion  suggested  that  the  freezing  of  supercooled 
drizzle  drops  increased  the  precipitation  efficiency  of  summertime  con¬ 
vective  clouds  in  the  Central  United  States  (Braham,  1964;  1968).  At 
the  time  that  these  ideas  were  introduced,  however,  quantitative  testing 
of  the  hypothesis  was  a  long  way  off.  Only  recently,  in  fact,  has  our 
knowledge  of  coalescence  and  riming  progressed  to  the  point  that  treat¬ 
ment  of  the  problem  could  be  attempted.  This  paper  performs  such  an 
evaluation,  examining  differences  in  particle  growth  by  riming  and 
coalescence,  and  seeing  whether  these  differences  would  be  maintained  in 
different  geographical  or  climatological  regimes. 

The  starting  point  for  the  current  discussion  is  the  familiar  growth 
equation  for  a  particle  of  mass  M,  -and  terni nal ••  veloci ty  V,  growing  by 
gravitational  sweepout  of  a  mcnodisperse  cloud  of  water  droplets  of 
mass  m  and  terminal  velocity  v: 

|r  =  A  U-v)  E  c  xc  ca  m  .  (1) 

whe^e  A  is  the  cross-sectional  area  of  tne  particle  normal  to  the  direc 
tier  of  fall,  >.  is  the  cloud  droplet  mixing  ratio,  p  is  the  air  den- 

L  3 

sity,  E  is  an  appropri ately  defined  collision  efficiency,  and  c  is  the 
fraction  cf  collisions  that  result  in  a  coalescence.  If  we  restrict 


ourselves  to  oblate  spheroidal  shapes,  the  cross-sectional  area  A  can 
be  expressed  as 


where  in  addition  to  the  terms  already  defined,  is  the  bulk  density 
of  the  collecting  particle  and  4  is  the  axis  ratio  of  the  oblate 
spheroid. * 


The  major  goal  of  this  study  is  to  examine  the  expected  growth  rates 
of  graupel  or  frozen  drops  as  compared  to  those  of  liquid  drops  of  the 
same  mass.  This  can  be  done  most  easily  by  normalizing  the  riming 
growth  rates  predicted  by  (1)  with  the  correspondi ng  coalescence  growth 
rate  for  an  unfrozen  drop  of  identical  mass.  This  results  in  a  quan¬ 
titative  estimate  of  the  relative  growth  rate,  to  be  represented  by  the 
symbol  y,  such  that 

v  _  dH  (graupel)  ,  dM  (liquid)  .  . 

dt  '  dt  ’ 

where  the  numerator  and  denominator  are  both  evaluated  using  (1).  The 
relative  growth  rate  defined  by  (2),  is  the  primary  quantity  that  will 
be  used  in  the  subsequent  discussions.  No  consideration  will  be  made  of 
ether  modes  of  growth,  such  as  condensation  or  aggregation.  In  many 
cases,  convective  precipitation  development  can  be  thought  of  as  a  two- 
stage  process  involving  production  of  embryos  that  are  large  enough  to 

*If  the  diameters  of  three  principal  axes  of  the  spheroid  are  Dj,  D^,  and 
D,  such  that  0}  =  03  and  D 3  <  D i ,  then  t  =  ^ 3 / D j . 


begin  to  grow  effectively  by  riming  or  coalescence  and  the  subsequent 
growth  of  these  embryos  to  precipitable  sizes.  This  study  focuses  on 
the  second  stage  of  that  process. 

2.  Formulation 

It  is  readily  apparent  that  there  are  three  critical  parameters  in 
(1)  that  need  to  be  evaluated:  the  terminal  velocity,  the  collision 
efficiency,  and  the  coalescence  efficiency.  Hi storical ly ,  each  of 
these  parameters  has  been  subject  to  considerable  uncertainty.  In 
each  case,  however,  recent  work  permits  us  to  make  more  quantatitive 
estimates  of  the  parameter  values  than  previously  possible. 

a .  Terminal  Velocities 

Terminal  velocities  for  water  drops  are  relatively  easy  to  estimate. 
Beard  (1976),  for  example,  gives  appropriate  expressions  as  functions  of 
drop  size,  ambient  temperature,  and  pressure.  Ice  particle  terminal 
velocities,  on  the  other  hand,  are  less  well  known  and  considerably  more 
complicated.  For  these  particles,  the  terminal  velocity  will  depend  on 
the  particle's  size,  shape,  bulk  density,  and  surface  roughness  as  well 
as  the  environmental  temperature  and  pressure. 

The  starting  point  for  the  graupel  terminal  velocity  estimates  will 
be  an  expression  for  the  fall  velocity  of  a  smooth  spherical  particle  of 
bulk  density  pj  falling  in  air  with  dynamic  viscosity  n  and  density  pa  . 

£ 

VS  =  TTi  I" 2  +  c?  ><  -  2  ( 1  +  C,  *>°-5l 


where 


M  end  D  are  the  particle  mass  and  diameter,  g  is  the  acceleration  of 
gravity,  and  Ci  and  C2  are  constants  equal  to  20.52  and  0.0901,  respec¬ 
tively.  The  dynamic  viscosity  of  the  air,  n,  is  exclusively  a  function 
of  temperature,  while  air  density,  pa,  depends  on  both  temperature  and 
pressure.  The  dimensionless  variable  is  often  termed  the  Best 
number.  Equation  (3)  is  ultimately  based  on  Abraham's  (1970)  expression 
for  the  drag  coefficient  of  a  sphere  as  a  function  of  Reynolds  number, 
as  discussed  by  Beard  (1980).  The  actual  terminal  velocity  of  the  ice 
particles  can  be  obtained  by  multiplying  V5  by  correction  factors  for 
shape  or  roughness: 


V  *  f  V$  .  (4) 

Since  a  particle's  drag  coefficient  is  inversely  proportional  to  the 
square  c*‘  its  terminal  velocity,  the  appropriate  correction  factor  will 
be  of  the  form 


£ 
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(^Ds/CV 


where  is  the  drag  coefficient  for  the  nonspherical  or  rough-surfaced 
particle  while  1S  the  drag  coefficient  of  a  smooth  spherical  par¬ 
ticle.  By  limiting  consideration  to  oblate  spheroids,  shape  effects  can 


be  described  as  an  additive  factor  of  the  form  Cq  =  Cqs  t  c,,  where  (,  is 
a  simple  function  of  the  axis  ratio  of  the  spheroid  (Beard,  1976).  In 


this  case,  £  =  0.3  ( 1 -  *  ) .  Roughness  effects  are  less  well  known,  but 
rr.ay  be  evaluated  by  com:  ring  graupel  drag  coefficient  estimates  based 
on  fall  velocity  and  mass  measurements  (e.g.,  Heymsfield,  1978)  with  the 
correspondi ng  coefficients  for  smooth  spheres.  Even  though  this  sort  of 
data  typically  includes  both  roughness  and  shape  effects,  there  is  a 
strong  suggestion  that  roughness  effects  may  be  expressed  as  a  multipli¬ 
cative  factor  of  the  form  Cg  =  cC^.  For  dry  graupel  particles  having 
Reynolds  numbers  >  100,  typical  surface  roughnesses  seem  to  correspond 
to  a  value  of  t  =  2.  For  Reynolds  numbers  <  100,  roughness  seems  to  be 
progressingly  less  important  and  ;  can  be  assumed  to  drop  linearly  with 
the  logarithm  of  the  Reynolds  number  from  2.0  at  Reynolds  number  100  to 
1.0  at  a  Reynolds  number  of  1.0.  If  surface  roughness  doubles  the  drag 
coefficient,  the  terminal  velocity  will  be  reduced  to  about  0.7  that  of 
a  water  drop  of  similar  mass,  which  is  exactly  the  reduction  in  fall 
speed  measured  by  Braham  (1964). 


it  may  be  interesting  to  note  that  the  "golf  ball"  effect,  in  which 
acded  roughness  in  the  form  of  dimpling  on  the  surface  of  the  ball 
reduces  drag,  is  not  applicable  to  frozen  raindrops  or  graupel.  In  the 
case  of  the  golf  call,  the  dimpling  serves  to  help  induce  turbulence  in 
the  boundary  layer  of  the  projectile  sooner  than  it  would  occur  if  the 
ball  were  smooth.  This  turbulence  delays  the  separation  of  the  boundary 
layer,  witn  a  correspond ng  reduction  in  drag.  This  phenomenon,  however 
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is  only  observed  at  very  high  Reynolds  numbers  (approximately  10^). 
Therefore,  if  it  is  applicable  to  all  atmospheric  particles,  it  is  only 
of  concern  for  large  hail  (Mason,  1971). 2 

b.  Col  1 i si  on  e  f  f i ci enci es 

A  cloud  droplet  being  overtaken  by  a  larger  particle  will  tend  to 
follow  the  airflow  as  it  diverts  around  the  particle.  This  means  that  a 
large  particle  will  actually  collide  with  only  a  fraction  of  the  cloud 
droplets  in  its  direct  path  of  fall.  This  fraction  is  termed  the  colli¬ 
sion  efficiency,  E.  In  general,  collision  efficiencies  between  falling 
particles  need  to  be  known  as  functions  of  the  size  shape,  bulk  density, 
and  surface  properties  of  the  interacting  particles,  and  as  functions  of 
the  ambient  air  temperature  and  pressure  as  well.  To  date,  most  studies 
of  collision  efficiencies  have  been  concerned  with  collisions  between 
liquid  drops  with  very  little  work  being  done  on  rime  growth.  The  major 
exceptions  to  this  generalization  are  the  wind  tunnel  study  of  Pflaum 
and  5ruppacher  (1979)  and  the  suggestion  by  Hall  (1980)  that  Beard  and 
Grover's  (1974)  algorithm  for  estimating  efficiencies  with  which  small 
raindrops  collide  with  cloud  droplets  or  large  aerosol  particles  could 
be  extended  to  include  eitner  liquid  or  solid  collectors.  In  their 
stucy,  Heard  and  Grover  followed  Langumir's  (1948)  approach  in  deve¬ 
loping  quantitative  expressions  for  the  collision  efficiencies  in  terms 
of  two  dimensionless  parameters:  the  Stokes  impaction  parameter  and  the 
Reynolcs  number  of  the  collecting  particle.  This  means  that  it  is  easy 

2 A  letter  in  June  1979  issue  cf  Scientific  American  has  an  interesting 
example  of  an  i nappropri ate  application  of  the  golf  ball  analogy  to 
small  particles. 
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to  obtain  collision  efficiency  estimates  for  almost  any  type  of  collec¬ 
tor  or  set  of  environmental  conditions  by  merely  making  appropriate 
modifications  to  these  two  parameters.  If  we  limit  our  consideration  to 
collectors  that  have  oblate  spheroidal  shapes,  the  Reynolds  number,  R, 
and  Stokes  impaction  parameter,  K,  can  be  expressed  as 


csc  S  d‘ 

9  o. 


where  in  addition  to  terms  already  defined,  C<.£  is  the  Cunningham  factor 
(see  Beard  and  Grover,  1974),  d  is  the  cloud  droplet  diameter,  and  p£  is 
the  density  of  liquid  water. 


This  approach,  however,  is  based  on  flow  around  smooth  spheres  and 
does  not  account  for  perturbations  in  the  flow  filed  caused  by  surface 
roughness  or  flow  through  the  particle.  Of  course,  surface  roughness 
still  plays  an  important  role  in  determining  the  collision  efficiency 
through  its  inclusion  in  the  terminal  velocity  estimates,  even  if  there 
is  no  attempt  to  consider  specific  disturbances  to  the  fluid  flow  intro¬ 
duced  by  c-e  roucn  surface.  In  particular,  increased  surface  roughness 
will  slow  tne  fall  of  the  graupel ,  reducing  its  characteristic  Reynolds 
numoer  and  this  in  turn  lowers  the  collision  efficiency. 


Hey ms  field  and  Pflaum  (1SS5)  recently  reexamined  Pflaum  and 
Pruppacher's  original  data  on  the  riming  of  frozen  drops  and  concluded 
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that  the  collision  efficiency  estimates  obtained  with  the  formulas  of 
Beard  and  Grover  generally  agreed  quite  well  with  the  wind  tunnel  esti¬ 
mates  of  mean  collision  efficiencies  during  graupel  growth.  While  the 
available  data  only  covered  a  portion  of  the  domain  of  interest,  this 
study  confirms  that  the  Beard  and  Grover  formulation  will  often  provide 
a  useful  first  approximation  for  the  collision  efficiencies  of  growing 
graupel.  Ultimately,  however,  it  is  likely  that  a  correction  term  will 
be  needed  to  adjust  these  calculated  efficiencies  for  factors  such  as 
the  shape,  bulk  density,  and  surface  roughness  of  the  graupel. 

For  simplicity,  collision  efficiencies  for  both  raindrops  and  graupel 
will  be  calculated  using  the  equations  of  Beard  and  Grover  (1974). 

While  other  sources  are  available  to  estimate  collision  efficiencies  for 
liquid  collectors,  the  Beard  and  Grover  formulation  tends  to  be  as  good 
as  any  of  the  alternatives  for  most  sizes  of  interest,  and  is  virtually 
unique  in  permitting  extrapolation  of  the  collision  efficiencies  to  dif¬ 
ferent  environmental  conditions.  In  addition,  it  is  clearly  advanga- 
gecus  to  use  the  same  basic  formulation  for  both  liquid  and  solid 
hydrometeors  to  avoid  any  bias  in  the  relative  growth  rates  that  might 
result  from  inconsistencies  in  the  collision  efficiency  calculations. 


In  calculating  graupel  growth  rates,  it  is  necessary  to  extend  the 
cal cu! at i cns  beyond  Reynolds  number  4Q0,  which  represents  the  upper 


limit  evlirtly  con 

si dered 

by  Beard  and 

Grove".  Beard  and 

Grover  did 

suggest  an  equation 

i  G  i  i  c  r 

cer  oart’cles, 

but  this  equation 

W3S 

merely 

intended  to  provide 

a  rough 

estimate  for 

large  collectors  using 

poten- 

tial  flow  approximations.  If  this  equation  is  arbitrarily  introduced 
for  all  R  >  400,  there  will  be  a  sharp  discontinuity  in  the  growth  esti¬ 
mates  at  R  =  400  since  the  potential  flow  calculations  result  in  signi¬ 
ficantly  increased  collision  efficiencies.  In  the  upcoming 
calculations,  the  transition  to  the  potential  flow  solution  is  smoothed 
by  simply  interpolating  between  the  two  formulations  over  the  Reynolds 
number  range  of  400  to  2400. 

c .  Coalescence  efficiencies 

The  coalescence  efficiency,  e,  represents  the  fraction  of  collisions 

* 

that  actually  result  in  coalescence.  Failure  to  coalesce  may  be  due  to 
the  drops  bouncing  apart  before  the  air  film  trapped  between  them  can 
drain,  or  by  the  rapid  breakup  of  a  large  drop  created  by  temporary 
coalescence  of  the  interacting  drops.  In  the  present  context,  the  first 
mechanism  is  the  one  of  primary  interest. 

While  it  has  long  been  known  that  colliding  drops  will  not  always 
coalesce,  it  has  been  difficult  to  develop  quantitative  expressions  for 
c  that  nave  general  appl i cabi I i ty .  Perhaps  the  best  available 
expressions  are  the  empirical  formulas  of  Beard  and  Ochs  (1984).  In 
using  these  equations,  however,  Ochs  and  Beard  (1985)  suggest  that  the 
calculated  efficiencies  be  constrained  tc  stay  in  the  ranee  of  50  to  100 
percent. 

In  analyzing  their  laboratory  data,  Beard  and  Ochs  (1984)  made  use  of 
a  modified  ’Weber  number.  This  dimensionless  parameter  involves  a 


characteristic  radius  of  deformation  and  a  representative  impact  speed 
(Schotland,  1950;  Foote,  1975).  In  Ochs  and  5eard  (19S4)  the  relevant 
deformation  radius  was  interpreted  as  the  small  droplet  radius  while  the 
characteristic  impact  speed  was  associated  with  the  terminal  velocity  of 
the  large  drop.  In  this  case  the  Weber  number  becomes 


where  p„  is  the  density  of  liquid  water,  d  is  the  diameter  of  the  small 

A, 

droplet,  V  is  the  terminal  velocity  of  the  collector  drop,  and  a  is  the 
surface  tension  of  liquid  water.  When  the  Weber  number  is  small,  the 
large  drop  does  not  deform  significantly  and  the  coalescence  efficiency 
is  near  unity.  As  the  Weber  number  increases,  the  possible  deformation 
of  the  large  drop  should  increase,  with  an  associated  reduction  'in  the 
coalescence  efficiency  (Ochs  and  Beard,  1984;  Beard  and  Ochs,  1984). 
This  analysis  supports  the  intuitive  conclusion  that  graupel  or  frozen 
drops ,  which  do  net  deform,  should  always  collect  the  supercooled  cloud 
dreps  t-ey  strike. 
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mined  a  variety  of  scenarios  involving  charged  drops  and  imposed 
electrical  fields,  there  is  no  generally  accepted  treatment  of  the 
problem.  As  a  result,  the  magnitude,  and  even  the  direction,  of 
possible  electrical  effects  remains  unresolved.  In  spite  of  these 
uncertainties,  there  is  usually  a  tendency  to  assume  that  electrical 
effects  may  be  expected  to  enhance  coalescence. 

3.  Results 

In  examining  the  relative  growth  rates  of  graupel  particles,  the 
range  of  parameters  considered  will  be  extended  well  beyond  the  con¬ 
ditions  Braham  observed  in  Missouri  clouds.  The  critical  variables 
involved  include  the  size,  shape,  bulk  density,  and  surface  roughness  of 
the  graupel  particle  as  well  as  the  diameter  of  the  cloud  droplets. 

Figures  1,  2,  and  3  illustrate  many  of  the  key  results  of  this  study, 
showing  the  relative  growth  rates  of  spherical  graupel  particles  of  two 
different  bulK  densities  as  functions  of  the  particle's  spherical  melted 
diameter,  0C,  where 


The  actual  diameter  of  tne  graupel  particles,  of  course,  would  be  simply 
1 " .  In  each  case,  the  curves  labeled  "smooth"  used  terminal 
velocities  for  tre  collecting  particle  based  on  drag  coefficients  for 
smooth  spheres  while  tne  curves  labeled  "rough"  used  the  surface  rough¬ 
ness  algorithm  for  d~y  graupel  discussed  in  section  2a.  Figure  I  was 
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calculated  for  a  cloud  droplet  diameter  of  12  pm,  while  the  figures  2 
and  3  are  based  on  cloud  droplet  diameters  of  16  and  20  pm,  respectively 
A  cloud  droplet  diameter  of  12  pm  may  be  appropriate  for  the  cold-based 
continental  clouds  of  the  High  Plains  (e.g..  Cooper  and  Lawson,  1984), 
while  the  larger  cloud  droplets  would  be  more  appropriate  for  warmer 
clouds  or  clouds  with  lower  droplet  concentrations. 

All  in  all,  these  results  support  Braham's  conclusion  that  the 
freezing  of  large  supercooled  drops  in  warm-based  clouds  may  signifi¬ 
cantly  accelerate  precipitation  development  (e.g.,  Fig.  3).  when  a 
supercooled  drizzle  drop  freezes,  there  should  be  an  immediate  jump  in 
the  accretional  growth  rate  associated  with  an  increase  in  coalescence 
efficiency  to  e  =  1.0.  Like  liquid  drops,  newly  frozen  drops  would  be 
expected  to  have  smooth  surfaces.  So  long  as  the  surface  remains 
smooth,  the  growth  rate  of  the  particle . wi 1 1  remain  high.  As  the  par¬ 
ticle  rimes,  however,  the  surface  will  quickly  roughen  (e.g.,  Keymsfield 
ana  Pflaum,  1985),  with  corresponding  reductions  in  the  terminal  veloci¬ 
ties,  efficiencies,  and  growth  rates  to  mare  modest  values.  In  this 
case,  crop  freezing  would  be  accompanied  by  a  sudden,  but  short,  spurt 
in  accretional  growth.  Even  after  this  spurt  nas  ended,  the  rime  growth 
"ates  may  still  be  si gni f icantly  larger  than  for  water  d-ops  cf  the  same 
~ass.  While  the  associated  growth  from  the  vapor  (sublimation)  would 
a1 so  be  greater  than  for  a  co- parable  unf-ozen  drop,  these  growth  rates 
are  usually  small  for  particles  of  this  size,  as  compared  to  growth 
rates  for  coalescence  or  riming,  and  may  oe  neglected. 


This  pattern  of  a  sudden  spurt  in  growth  following  freezing  would 
also  be  expected  in  clouds  composed  of  smaller  droplets  (e.g.,  figs.  1 
and  2).  The  main  difference  would  be  that  the  magnitude  of  the  growth 
spurt  will  be  smaller,  and  the  subsequent  growth  rates  may  often  slow  to 
less  than  that  of  a  liquid  drop  of  comparable  mass. 

One  of  the  interesting  features  of  this  analysis  is  the  relative 
insensitivity  of  the  relative  growth  rates  to  changes  in  the  bulk  den¬ 
sity  of  the  graupel.  This  result  is  examined  in  more  detail  in  figure 
4,  illustrating  the  relative  growth  rates  as  functions  of  the  graupel 
bulk  density.  In  each  plot,  separate  curves  are  drawn  for  D0  =  0.5  mm 
and  D0  =  2.0  mm.  In  these  examples,  only  the  surface  roughness 
algorithm  for  dry  graupel  was  used. 

For  a  cloud  droplet  diameter  of  12  um,  the  growth  curves  are  relati¬ 
vely  flat  with  a  hint  of  a  maximum  in  the  bulk  density  range  of  0.5  to 
0.6  g  cm~3  before  dropping  off  rather  sharply  as  the  bulk  density  falls 
below  0.1.  For  slightly  larger  cloud  droplets  (d  =  16  um),  the  relative 
growth  rates  are  larger  for  all  bulk  densities,  with  maximum  values 
associated  with  bulk  densities  in  the  range  of  0.2  to  0.3  g  cm"3.  For 
d  =  20  um,  the  relative  growth  rates  increase  still  further  with  the 
maximum  values  corresponding  to  bulk  densities  near  0.1  g  cm"^.  In  the 
present  model,  the  treatment  of  very  low  bulk  density  particles  (<  0.1  g 

_3 

cm  )  is  suspect  since  no  airflow  through  the  particle  is  permitted. 
While  it  is  not  necessarily  clear  what  the  net  effect  of  such  flows 
would  be,  it  seems  that  at  the  very  least  this  would  tend  to  increase 
the  particle  terminal  velocities  and  result  in  increased  collision 
efficiencies  and  accretional  growth  rates. 


Up  to  this  point,  only  spherical  graupel  particles  have  been  con¬ 
sidered.  Figure  5  extends  the  calculations  to  oblate  sphenoidal  shapes 
with  axis  ratios  down  to  0.3.  As  in  figure  4,  these  examples  all  use 
the  surface  roughness  algorithm  for  dry  graupel.  And  again,  as  in 
figure  4,  separate  curves  are  drawn  for  D0  =  0.5  mm  and  D0  =  2.0  tim.  In 
this  case  the  maximum  graupel  diameter  would  be  simply  Dq  [pt/(*p-)]^ 
All  the  plots  in  this  figure  were  based  on  a  cloud  droplet  diameter  of 
16  pm,  with  separate  panels  correspond!' ng  to  bulk  densities  of  0.1,  0.3, 
and  0.7  g  cm" 3.  While  there  is  an  interesting  crossover  in  the  curves 
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in  each  panel,  the  principal  result  is  that  the  relative  growth  rates  do 
not  show  great  variation  with  shape.  For  the  warm-based  clouds, 
however,  there  is  a  distinct  tendency  for  increased  growth  rates  as  the 
axis  ratio  gets  smaller. 

Since  the  collision  efficiency  estimates  are  based  on  spherical 
collectors,  the  predicted  collision  efficiencies  are  .likely  to  be 
underestimated  for  low  axis  ratio  oblate  spheroids.  Hall  (1980)  has 
developed  an  alternative  set  of  equations  to  estimate  collision  effi¬ 
ciencies  for  "plate-like"  particles  which  can  be  used  to  illustrate  the 
-arce  cf  uncertainty  present  in  calculated  growth  rates  for  4 <<  1. 
Hall's  collision  efficiency  formulation  parallels  that  of  Beard  and 
Grover  (1974),  but  is  based  on  the  results  of  Fitter  and  Pruppacher 
(1974)  and  Pitter  (1977)  for  oblate  spheroids  of  axis  ratio  0.05.  As  in 
the  work  of  Beard  and  Grover,  the  collecting  particles  are  assumed  to 
have  smooth  surfaces  that  are  impermeable  to  airflow.  Relative  growth 
rates  using  Hall's  collision  efficiencies  are  plotted  as  dashed  lines  in 
figure  5. 


As  should  be  expected,  the  relative  growth  rates  are  significantly 
larger  using  Hall's  collision  efficiencies.  In  addition,  these  results 
do  not  evidence  the  low  bulk  density  drop  off  in  growth  rates  predicted 
for  spherical  collectors  (e.g.,  Fig.  4).  Although  the  current  model  is 
not  well  suited  to  evaluate  the  riming  of  snowflakes  and  aggregates,  the 
results  suggest  that  these  low  density  particles,  which  often  have  axis 
ratios  in  the  range  of  0.3  to  0.5,  may  be  quite  effective  growth  centers 
for  accretion. 3 

4.  Discussion 

This  analysis  of  the  relative  growth  rates  of  graupel  and  supercooled 
raindrops  tends  to  confirm  Braham's  conclusion  that  graupel  particles  in 
warm-based  clouds  may  be  expected  to  grow  faster  by  riming  than  unfrozen 
drops  of  the  same  mass  grow  through  coalescence.  The  Missouri  clouds 
that  Braham  was  studying,  however,  seem  particularly  favorable  to 
graupel  growth  and  his  conclusions  should  not  be  generalized  to  clouds 
in  other  geographical  areas,  such  as  the  High  Plains.  In  fact,  if  the 
supercooled  cloud  droplets  are  small  enough,  graupel  particles  can  often 
grow  significantly  slower*  than  waterdrops  of  comparable  mass.  These 
microphysical  growth  rates,  of  course,  are  only  one  aspect  of  precipita¬ 
tion  efficiency.  Since  the  changes  in  accretional  growth  rates  asso¬ 
ciated  with  phase  changes  can  be  ratner  dramatic,  however,  it  seems 


^Recent  studies  oy  Rasmussen  an.i  Lew  (1935)  also  suggest  that  aggregates 
may  rime  more  effectively  than  previously  thought. 


likely  that  this  could  well  influence  the  probability  that  a  cloud  deve¬ 
lops  precipi t at i on  or  the  amount  of  precipitation  produced. 

In  examining  the  relative  growth  rates  defined  by  (1)  and  (2),  it  is 
apparent  that  graupel  particles  are  increasingly  favored  over  liquid 
drops  as  the  collecting  particles  get  larger.  Increasing  the  cloud 
droplet  diameter,  d,  can  also  produce  significant  increases  in  the 
relative  growth  rates  and  is  probably  the  critical  factor  in  controlling 
the  geographical  variability  of  the  present  results.  In  general,  warm 
cloud  bases  and  low  droplet  concentrations  promote  large  cloud  droplets 
while  cold  cloud  bases  and  high  droplet  concentrations  tend  to  reduce 
the  mean  droplet  diameter.  The  sensitivity  to  droplet  size  is  linked  to 
the  coalescence  efficiency,  which  is  one  of  the  two  critical  terms  that 
determine  the  difference  between  the  growth  of  liquid  and  solid  collec¬ 
tors.  The  importance  of  this  term  to  the  overall  result  also  suggests 
that  the  charged  particles  and  electric  fields  of  a  mature  thunderstorm 
may  often  affect  the  mi crophysi cal  efficiency  of  precipitation  develop¬ 
ment.  For  example,  if  the  charged  environment  in  a  thunderstorm  were 
found  to  increase  coalescence  efficiencies  to  values  near  1.0,  raindrop 
coalescence  growth  might  be  enhanced  by  as  much  as  a  factor  of  two. 

It  is  important  to  note  that  the  increased  growth  expected  for  a 
newly  frozen  raindrop  as  c  +  1  is  seldom  realized  for  very  long.  In 
general,  this  increase  will  progressively  be  offset  oy  the  increasing 
surface  roughness  as  the  particle  rimes,  slowing  its  fall  and  lowering 
the  collision  efficiency.  The  end  result  is  thus  a  balance  between  the 
effects  of  different  coalescence  efficiencies  and  surface  roughness  fac- 
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tors  for  liquid  drops  and  graupel.  Other  factors,  such  as  particle 
shape  and  bulk  density,  also  play  a  role  in  the  microphysical  growth 
rates,  but  are  generally  less  important. 
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The  insensitivity  to  changes  in  particle  bulk  densities  is  par¬ 
ticularly  interesting  to  examine.  From  (3),  it  is  apparent  that  for 
spherical  particles  of  the  same  mass  and  roughness,  the  fall  velocity  is 
essentially  proportional  to  the  inverse  of  the  particle  diameter.  The 
area  swept  out,  on  the  other  hand,  is  directly  proportional  to  the 
square  of  the  diameter.  Decreasing  the  bulk  density  thus  increases  the 
volume  swept  out  by  the  graupel  particle,  in  an  amount  proportional  to 
the  particle's  increase  in  diameter.  The  lower  fall  velocity,  however, 
reduces  the  collision  efficiency  and  almost  totally  offsets  the 
increased  sweep-out  volume.  This  balance  between  increasing  the  sweep- 
out  volume  and  decreasing  the  collision  efficiency  is  also  the  key  fac¬ 
tor  in  explaining  the  lack  of  sensitivity  to  changes  in  shape. 

The  most  direct  application  of  this  study  is  to  compare  particle 
growth  before  and  after  a  change  of  phase.  This  could  involve  the 
freezing  of  a  supercooled  drizzle  drop  or  melting  of  a  falling  graupel 
particle  before  it  is  recycled  into  another  updraft.  In  addition,  these 
analyses  can  be  applied  to  wet  growth,  where  a  graupel  particle  collects 
supercooled  water  so  quickly  that  the  latent  heat  released  cannot  be 
dissipated  rapidly  enough  to  freeze  all  the  water  that  is  accreted.  In 
most  cases,  wet  growth  should  generate  a  smooth  surface  on  the  growing 
graupel  with  an  associated  jump  in  growth  rate.  Examination  of  figures 
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1  to  3,  in  fact,  suggests  that  the  onset  of  wet  growth  may  often  mean  a 
doubling  in  the  growth  rate.  This  means  that  the  transition  back  to  dry 
growth  will  often  take  place  at  substantially  lower  water  contents  than 
required  to  initiate  wet  growth  in  the  first  place.  Since  such  tran¬ 
sitions  depend  on  the  prior  history  of  the  particle,  they  represent  a 
good  example  of  a  hysteresis  loop.  At  some  point,  as  the  water  coat 
becomes  thicker,  the  surface  of  the  wet  graupel  should  start  to  deform 
during  collisions  with  the  result  that  the  coalescence  efficiency  may 
begin  to  approach  that  of  an  unfrozen  waterdrop. 

This  study  can  also  be  used  to  compare  the  relative  growth  rates  of 
different  kinds  of  precipitation  embryos.  It  is  necessary  to  keep  in 
mind,  of  course,  that  all  types  of  embryos  are  not  equally  likely  to  be 
found  in  any  given  cloud  and  thus  the  comparisons  are  largely  academic. 
In  addition,  it  is  not  imperative  that  relative  growth  estimates  of  this 
sort  be  made  on  the  basis  of  particles  of  the  same  mass.  Heymsfield 
(1932),  for  example,  has  performed  a  similar  comparison  for  particles  of 
the  same  diameter.  The  growth  rate  is  so  dependent  on  the  particle 
~>ass,  however,  that  the  results  of  such  comparisons  simply  show  that 
particles  of  high  bulk  density  grow  faster  than  comparable  diameter  par¬ 
ticles  of  low  bu!<  density.  The  interesting  feature  of  the  present  ana¬ 
lysis  is  that  the  bulk  density  depenaence  is  reduced  to  a  second  order 
effect  when  tne  growth  rate  comparisons  are  made  on  the  basis  of  par¬ 
ticles  of  the  same  mass.  In  fact,  the  lower  bulk  density  graupel  forms 
were  often  found  to  be  slightly  more  effective  as  precipitation  embryos 
than  higher  bulk  density  graupel  or  frozen  drops. 


While  the  overall  pattern  of  the  results  should  be  fairly  accurate, 
it  is  important  to  remember  that  both  of  the  critical  parameters  that 
control  the  magnitude  of  the  relative  growth  rates  are  still  subject  to 
considerable  uncertainty.  Our  quantitative  knowledge  of  coalescence 
efficiencies,  for  example,  is  still  in  its  infancy  and  almost  no  quan¬ 
titative  work  has  been  done  on  the  effects  of  surface  roughness.  If 
anything,  however,  these  lingering  uncertai nti es  only  highlight  Braham's 
insight  in  his  pioneering  analysis  of  Missouri  clouds. 
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Figure  Captions 


Figure  1.  Relative  growth  rates  of  graupei  particles  with  bulk  densities 

of  0.3  and  0.7  g  ca  ^  as  functions  of  celted  drop  diameter  ( ) . 

In  this  case,  the  results  are  based  on  a  cloud  droplet  diameter 
of  12  pa. 

Figure  2.  Relative  growth  rates  of  graupei  particles  with  bulk  densities 

_3 

of  0.3  and  0.7  g  ca  as  functions  of  celted  drop  diameter  ( Dq ) . 

In  this  case,  the  results  are  based  on  a  cloud  droplet  diameter 
or  16  pm. 

Figure  3.  Relative  growth  rates  of  graupei  particles  with  bulk  densities 

-3 

of  0.3  and  0.7  g  cm  as  functions  of  celted  drop  diameter  (D^), 

In  this  case,  the  results  are  based  on  a  cloud  droplet  diameter 
of  20  pm.  , 

Figure  4.  Relative  growth  rates  of  rough  surface  graupei  particles  of  melted 
diameters  Q.5  and  2.0  mm  as  functions  cf  particle  bulk  density. 
Separate  plots  are  shown  fer  cloud  droplet  diameters  of  12,  16,  and 
20  pm. 


Figure  5.  Relative  growth  rates  of  rough  surface:  oblate  spheroidal  graupei 

particles  of  celted  diameters  0.3  and  2.0  cm  as  functions  the  axis 

ratio  of  the  particle.  The  dashed  l*'r.es  at  the  left  of  each  panel 

are  based  on  collision  efficiencies  for  plate-like  particles  (see 

text  for  discussion).  Separate  plots  are  shewn  for  bulk  densities 

-3 

of  0.1,  0.3,  and  0.7  g  cm  .  In  each  case,  the  results  are  based 
on  a  cloud  droplet  diameter  of  16  „c. 
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Figure  2.  Relative  growth  rates  of  graupel  particles  with  bulk  densities 


of  0.3  and  0.7  g  cm  as  functions  of  melted  crop  diameter  (D  ) 
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In  this  case,  the  results  are  based  on  a  cloud  droplet  diameter 


or  16  pm. 
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Figure  5.  Relative  growth  rates  of  rough  surfaced  oblate  spheroidal  graupel 
particles  of  celted  diameters  0.5  sr.d  2.0  na  as  functions  the  axis 
ratio  of  the  particle.  The  dashed  lines  at  the  left  of  each  panel 
are  based  on  collision  efficiencies  for  plate-like  particles  (see 
text  for  discussion).  Separate  plots  are  shovn  for  bulk  densities 

_3 

of  0.1,  0.3,  and  0.7  g  cm  .  In  each  case,  the  results  are  based 
on  a  cloud  droplet  diameter  of  16  pm. 
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